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ABSTRACT 


Process control models can be developed by using both conventional and AI 
approaches. The conventional approaches include regression models and process 
kinetic models. The artificial intelligence (AI) approaches are based on artificial 
neural nets (ANN), genetic algorithm (GA) and fuzzy rule based expert systems 
(FREES). Plant data on hot metal desulphurization, carried out by injecting calcium 
carbide in 400 ton torpedo ladles, is analyzed to test and tune different types and 
combinations of models and then evaluate their relative performance. While the 
control models based on process fundamentals provide a fillip to the new ideas for 
technological developments and improvements, the combination of conventional and 
AI approaches may be a better option for process control on the shop floor. It is 
advisable to first develop and test all models and then decide about the best strategy 
of using them. 

The phenomenon of nitrogen and oxygen pickup by the liquid steel from the 
atmosphere, during tapping into a ladle, is described first. The results of the 
absorption model are fed into the model for inclusion control. The kinetic model for 
inclusion control (prediction of inclusion size, size distribution and amount) is based 
on mixed mass transfer of dissolved aluminum and dissolved oxygen in the bath, 
collision mechanisms (Brownian motion. Stake’s collision, and turbulent collision), 
and simultaneous assimilation of inclusions by the top slag. The model also predicts 
the amount of aluminum needed to achieve a desired level of dissolved oxygen in the 
final steel. In the user-friendly model developed in this work, an output window has 
been created to view all the predictions in the form of graphic displays. Computation 
time has been optimized by applying genetic algorithm. 



The model has been tested for the case of 300 ton ladles at Bokaro Steel Plant. The 


SEM observation of the inclusions in the samples collected at Bokaro Steel Plant 
show that most inclusions, in absence of bottom stirring of ladle are of size 15 
microns or smaller. Model predictions agree well with this result. Further, according 
to the models developed in this work, only a maximum of 8% of the inclusions may 
get removed during a holding/transfer time of 45 minutes. Gas stirring in the ladle at 
the rate of 2 liters/ton/min, soon after tapping, would greatly help in the floatation of 
inclusions by more than 50% during the time of transfer to ladle treatment station. 
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jcinetics. The present work focuses on comparing the efficacy of different control 
models that can be developed for a given plant situation and the strategy that should 
be adopted in a plant environment. 

We shall first describe, briefly, various models, namely an empirical (regression) 
model (section 2. 1 .1), a simple kinetic model based on second order chemical reaction 
control (section 2.1.2), a sophisticated kinetic model based on the assumption ol'mass 
transfer (first order) control (section 2. 1 .3), ANN models (section 2. 1 .4), fuzzy rule 
based expert systems (FREES) model (section 2.1.5), and GA optimized model 
(section 2.1.6) to control/predict desulphurization of hot metal. Then we shall 
compare their efficacy in terms of the success to predict the final sulphur content of 
the hot metal. For the sake of a just comparison, the same set of data (section 2.2) is 
used to validate different models because the main point here is the importance of 
comparative performance of the models for a given plant situation. 

The liquid steel, at the end of heat in a basic oxygen furnace (EOF), may contain 400- 
900 ppm dissolved oxygen, 2-4 ppm dissolved hydrogen, 40-60 ppm dissolved 
nitrogen, and the temperature of metal varies in the range 1873-1973K. During 
tapping from EOF, the falling stream of liquid metal as well as the turbulent surface 
of liquid steel in the ladle promotes absorption of additional oxygen, hydrogen and 
nitrogen from the ambient atmosphere into the melt. After the tapping operation is 
over, deoxidizers like coke, Fe-Si, Fe-Mn, A1 may be added during the tapping 
operation and at the steel refining unit (SRU). For example, the steel may be subjected 
to argon stirring at argon rinsing station (ARS), or steel may be reheated to a desired 
temperature in the ladle furnace (LF) and thereafter subjected to further deoxidation 
treatment by aluminum and calcium. During the time of tapping and the time of 
transfer of ladle to SRU, the deoxidation reactions, growth of deoxidation products in 


2 



size due to mass transfer and collisions, and assimilation of deoxidation products into 
top slag take place continuously. The liquid steel in the BOF, during tapping and 
transfer to SRU, continuously looses heat by radiation and convection. The slag lying 
on top of metal can react with metal and with the lining of the ladle. For a better 
process control, it is necessary to understand the thermodynamics and kinetics 
occurring at various stages of steel processing. 

The product of deoxidation by aluminum, when aluminum is added to liquid steel 
during the tapping operation, is solid alumina. The presence of solid alumina in liquid 
steel above a certain level is undesirable, because alumina can clog the nozzle in 
tundish. Retention of solid alumina, as inclusions in the solidified product, is also 
harmful because it decreases cleanliness of the steel and affects the final mechanical 
properties of the product. It is necessary to control both the amount of size range of 
inclusions in liquid steel. During their residence time in liquid, the oxides in solid or 
liquid state, can grow in size by several mechanisms ' ’■* such as growth by mass 
transfer, Oswald ripening, and collision. The bigger inclusions (10-100 microns) can 
float up whilst the smaller inclusions (<10 microns) may get recirculated with the 
convection streams present in the liquid metal. As a thumb rule, the inclusions of size 
greater than 50 microns should be completely avoided in Hat products because they 
act as the sites for crack initiation. In thin sheets, as a precautionary measure, the total 

I 

amount of inclusions should be less than 30 ppm, 

The present study focuses on killed steels made in 300 ton BOFs at Bokaro steel 
plant. In the operating practice followed at Bokaro Steel Plant, the liquid steel is 
deoxidized in the ladle during tapping by adding deoxidizers such as carbon in the 
form of coke, ferro-silicon, ferro-manganese and aluminum. The amount of 
deoxidizer(.s) to be added is decided by the aim dissolved oxygen content of liquid 


3 



steel after tapping from BOF. Coke is added in the very initial stages ui Lapping su ac 
to remove part of the dissolved oxygen as CO gas. While the CO gas escapes from 
liquid metal almost instantaneously, the oxide products (Si02, MnO, AI2O3) may take 
finite time to float up. Chapter 3 describes the absorption of nitrogen and oxygen 
gases during tapping, and thermodynamic aspects of deoxidation by aluminum are 
discussed first, followed by the kinetic model of growth of alumina inclusions. 

It is reported in literature that the growth of inclusions may take place by several 
mechanisms like mass transfer. Brownian motion, turbulent collision and Stokes . 
collision. In this work, different models of inclusion growth have been developed and 
subjected to sensitivity analysis. A model, with a user friendly interface, is presented 
for the situation at Bokaro Steel Plant with the help of which it is possible to predict; 

• Pickup of nitrogen and oxygen from ambient atmosphere during the tapping 
operation (section 3.5.1), 

• Initial amount of deoxidant to be added for a given composition and 
temperature of steel (section 3.5.2)i, 

• Size distribution and amount of inclusion present in liquid steel as a function 
of time (section 3.5.4 and 3.5.5). 


4 


S 


i 



Chapter - 2 

2. Strategies for Development of Process Control 
Models for Hot Metal Desulphurization: 
Conventional and AI techniques 

2.1 Different models for prediction of final sulphur content of the bath 

There have been several studies^^'^'^’^^^ in the last decade to develop control/prediction 
models for hot metal desulphurization on the basis of metallurgical kinetics 
statistics (regression models) artificial neural nets (ANN) genetic algorithms 
(GA)''^'^^\ and fuzzy rule based expert systems (FRBES) The set of computational 

techniques comprising of ANN, GA and fuzzy logic (including their combinations 
like Fuzzy-GA, GA-ANN, Fuzzy-ANN, etc.) are now known as ‘soft computing 
techniques’. The different models are described below 

2.1.1 Statistical analysis of data and development of regression models 

The first step in model development should be to study the mean, variance, normal 
distribution, variance-covariance and correlation matrices, eigen values and eigen 
vectors followed by principal components analysis so as to decide the possibility of 
dimensionality reduction and check the presence of clusters in the data. Each column 
of data (for 229 heats, typical data set given in Table 2.1) was checked for normal 



distribution. By observing the magnitude of the off-diagonal correlations in the 
correlation matrix (given below) it was concluded that the variables (1) treatment 
time, (2) mass of hot metal, (3) initial sulphur content and (4) CAD injection rate, 
assumed to be independent, were only weakly correlated.. 


Correlation Matrix 


Treatment 

Mass of 

Initial S 

CaD injection 

Finals 

time 

hot metal 

content 

rate 

content 

1.0000 

0.2304 

0.4160 

-.2021 

-.6712 

0.2304 

1.0000 

0.0008 

-.0797 

-.0519 

0.4160 

0.0008 

1.0000 

-.0207 

0.2052 

-.2021 

-.0797 

-.0207 

, 1.0000 

-.0685 

-.6712 

-.0519 

0.2052 

-.0685 

1.0000 


On analyzing the eigen values, eigen vectors followed by principal components 
analysis, it was formd that there were no clusters in the data. The dimensionality of 
the problem could not be reduced because the scree curve (Fig. 2.1) did not show any 
rocky debris or flat portion at tail end. Multiple linear regression analysis was, 
therefore, carried out with final sulfur content as the dependent variable and the other 
four (described above) as independent variables. The standard error of estimate was 
0.0018 and multiple correlation coefficient (R) was . 0.91 (R squared = 0.8229, 
corrected R squared =0.8202). 




0 1 2 3 4 5 6 

f 

eigenvalue number 

Fig. 2. 1 : Scree-cruve: eigen value vs. eigen value number; variable 1 : Treatment time, 
2. Mass of metal, 3. Initial Sulphur content, 4. CaD injection rate, 5. Final Sulphur 


content. 




Several modifications can be made in the choice of input (independent) variables to 
be used in a regression model. In order to account for any possible non-linear 
correlation, the independent variables can be suitably modified. In the case of 
desulphurization of hot metal, for example, the independent variables can be suitably 
modified and selected on the basis of process fundamentals. As will be discussed 
later, according to kinetic fundamentals, the injection rate is expected to be dependent 


on initial and final sulphur in the form In 




. A linear regression model, called a 


ballistic model, was developed at Sollac, France, according to which the amount of 


powder injected, Qcao, was calculated from 
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^ /■" 


QcaD = Xi + Wm X In 






X {X2 + X3 X Sj + (X4 + X5 X Ir) X Sf + Xs X Ir} (2.1) 


Where, Xi X^: Constants, Wn,: Hot metal weight, Si : Initial sulfur content, Sf: 

Final sulfur content, Ir: Injection flow rate, Qcan: Mass of calcium carbide needed. 


2.1.2 Model based on second order chemical reaction 

In this model the desulphurization is assumed a priori to be controlled by second order 
chemical reaction and the governing equation is: / 

— = -KS,C (2.2) 

dt 

Where, C is the concentration of desulphurizer at time t, St is concentration of sulphur 



at time t and K is the chemical rate constant for desulphurizer. It is further assumed 
that the residence time (t) of any desulphurizer particle is short, relative to the total 
injection time, and during this injection time the hot metal sulphur content remains 
approximately constant. Upon integration of equation (2.2) final expression is 
obtained in terms of the injection rate (kg of powder/metric ton of hot metal) as 



(Sj-S,) 


MPf + 


MPf 


KT(Si-S,) 




(2.3) 


Where, M - - — —x 0.1, Ms and McaC 2 are molecular mass of sulfur and 

^CaC2 

desulphurizer respectively, P is the component purity, F is the blend fraction, Sj: 
initial sulfur content, Sf: sulfur content in the hot metal after time t. From the plot of 
actual versus predicted sulphur, on using the value of K as 8.33 and x = 1 s, the value 
of is obtained as 0.80. In this kinetic model no separate treatment is incorporated to 
account for the specific contributions of transitory reaction (while the particle is 
traveling in the metal) and the permanent contact reaction (when the particle has 
joined the top slag reactions). 


2.1.3 Model based on mass transfer (first order), including the transitory and 
permanent reactions 

A detailed kinetic model, considering the specific contributions of the transitory and 
the permanent contact reactions, assuming mass transfer as rate controlling step, is 
described in detail elsewhere This model was further improved by Sheshadri et al 
According to this model the desulphurization process is divided into three parts: 



contribution of individual powder particles while they are rising independently in the 
liquid metal, Rp, contribution of the powder particles which are captured in the rising 
carrier gas bubbles, Rt, and the contribution of particles when they join the top slag, 
Rts- The overall rate equation is 


d[%S] 

Mg dt 


-Rts +Rp +Rb 


(2.4) 
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2.1.4 ANN Model 

A detailed discussion of the application of neural nets to the desulphurization problem 
is explained elsewhere The artificial neural nets, similar to the regression 
models, can be directly applied to raw data without any prior knowledge of the 
process fundamentals. The sigmoidal function used for training ANN is able to 
automatically account for any nonlinearity present between the independent and 
dependent variables, and there is no need to hunt for a specific relationship similar to 
equation (2.2). However, better results can be obtained if the selection of the input 

and output variables is guided by prior process knowledge. 



The list of various symbols and their values used in equation (2.4) is given below. 


Symbol 


Value 

As-m 

Area of slag-metal interface (m'^) 

40 

db 

Average bubble diameter (m) 

Variable 

dp 

Average particle diameter (m) 

15 X 10’" 

f 

Fraction of particles entrapped in the bubbles 

0.3 

hi 

Depth of injection of CAD (m) 

1.8 

Ir 

Rate of injection of the desulfurizing agent (kg/s) 

0.45-1.5 

kb 

Sulfur partition coefficient between bubble 

' 

2.36x10'^ 

kf 

Sulfur transfer coefficient through the metal slag 

8.70 X 10"^ 

kp 

Sulfur transfer coefficient at particle melt boundary layer 

2.36 X 10'^ 

Lg'’ 

Sulfur partition coefficient between particles and melt, 

respectively. 

2000 

Ls® 

Sulfur partition coefficient between particles entrapped 

in the bubble and melt. 

2000 

L 

Sulfur partition coefficient between top slag and melt 

2000 

McaCz 

Molecular mass of CaCg (g/mol) 

64 

Ml 

Weight of metal (Ton) 

280 -370 

Ms 

Molecular mass of sulfur (g/mol) 

32 

Mts(t=o) 

Mass of carry over slag (kg) 

5 

Po 

Average pressure in side the bubble (atm) 

1.0 

Pr 

Average reference pressure (atm) 

Variable 

Q 

Gas flow rate (Nm^/s) 

0.033 

[%S] 

Instantaneous concentration of sulfur in the metal bulk 

Variable 




(%wt) 





Sf 

Final Sulphur content (%wt) 

0.002-0.02 

[%S]i 

Instantaneous concentration of sulfur at top slag- metal 

interface (metal side) (%wt) 

Variable 

[%S]o, 

Si 

Initial sulfur concentration in the melt (%wt) 

0.0015-0.07 

t 

Total time for desulphurization treatment (s). 

500-3400 

trp, T 

Particle residence time (s) 

Variable 

trb 

Bubble residence time (s) 

Variable 

PL 

Specific weight of melt (kg/m^) 

7000 

Pp 

Specific weight of particle (kg/m^) 

1740 




In the ANN model, a multidimensional surface is represented by a complex 
combination of terms, for example exponential equations when a sigmoidal function 
used for training. It is necessary that enough data is available to train the net for each 
part of the multidimensional surface, other wise in the regions where data is scarce the 
predictions can become very erroneous. This does not usually happen with other 
models described above which are based on fundamental considerations of process 
thermodynamics and kinetics. For this reason, it is advised to use ANN models for 
interpolation and not for extrapolation. 

The popular back propagation method was used in the present work for the training by 
employing a sigmoidal function, with momentum factor a = 0.5, training rate P = 0.5, 
and logistic T = 1.0. The training set consisted of 91 data and test set contained 79 
data. The training was terminated when the learning curve became nearly flat 
(visually); maximum iterations needed were less than 5000. Several net architectures 
were tried. The values obtained from the plot of actual versus predicted sulphur for 
some of nets, which yielded good results, were: for three layer (4-5-1 net) 0.87, three 
layer (4-6-1 net) 0.81, four layer (4-6-3-1 net) 0.80, and four layer (4-5-3-1) net 0.86, 
respectively. 

It is our experience that for most iron and steel making processes (viz. blast furnace, 
basic oxygen furnace, RH degasser, etc) if there are, say, n input variables and one out 
put variable, then a simple three layer net of n+1, n-1, n+2 or n-2 middle layer 
neurons usually gives the best results. The minimum number of training data needed 
may vary from 20n to 40n, but it should be properly distributed over the entire 


multidimensional domain. 



2.1.5 Fuzzy rule based expert system (FREES) model 


In principle, the fuzzy rule based expert system models are normally used when the 
fundamental knowledge of the process is poor, the relationship between the variables 
is complex and some of the process variations are stochastic. With regard to the 
desulphurization treatment, the ambiguity or imprecision in information is usually 
present in the measurement of input and output data coupled with stochastic variations 
in particle shape, size and composition, and fluid flow patterns with in the metal. For 
example, the concentration of sulphur measured at one location is not representative 
of the other parts of the same system. The concentration at the top of bath is different 
from that at the bottom due to improper mixing. Thus, even though concentration has 
been measured to an accuracy of ±1%, the system is fuzzy (or dizzy) with regard to 
the concentraton in different parts of the same system. As the vessel size changes, viz 
due to erosion of refractory, the flow patterns also change and process performance 
drifts. In the FREES approach, an attempt is made to divide the multidimensional 
surface of the input and out put variables into a number of regions demarcated by 
process and plant specific fuzzy decision rules so as to capture the process behavior. 
In our case, the output variable is final sulphur content of hot metal and the input 
variables are (1) treatment time, (2) mass of hot metal, (3) initial sulphur content of 
hot metal and (4) CAD injection rate. A lot of practical experience or expertise is 
required in fra min g the rules, and hence they can vary from person to person and from 
one plant to another. The obvious disadvantage is that a slight change in the operating 
conditions may call for the review of the whole set of rules framed to predict the 
output. On the other hand, an advantage of the fuzzy systems is that they can produce 

a nonlinear functional mapping just as the algorithms do. 



In the FREES model developed in the present work, treatment time is hizzified into 
five classes, CAD injection rate into four classes, hot metal weight into five classes, 
initial sulphur content into five classes, and final sulphur content into five classes. The 
linguistic labels used are: Very low (VL), Low (L), Medium (M), High (H) and Very 
high (VH). S3mimetrical, normal, triangular membership functions with 50% overlap 
are assigned for each fuzzy subset. The symbol p is used to represent the membership 
value. A typical example of the mathematical description of membership function of 
treatment time and its validity range is given in Table 2.2. The eleven fuzzy mles, 
given below, were arrived at after hundreds of trials with respect to A: treatment time, 
B: hot metal mass, C: initial mass % of sulphur, D: CAD flow rate, and E: final mass 
% sulphur in hot metal. 

1. IFAisVLandBisMandCisMandDisVHTHENEisVH. 

2. IF A is VH and B is L and C is M and D is VH THEN E is VL. 

3. IF A is VH and B is M and C is L and D is VH THEN E is VL. 

4. IF A is VL and B is M and C is VL and D is H THEN E is L. 

5. IF A is H and B is VH and C is M and D is VH THEN E is L. 

6. IF A is H and B is H and C is M and D is H THEN E is M. 

7. IF A is VL and B is M and C is L and D is H THEN E is H. 

8. IF A is M and B is M and C is M and D is VH THEN E is H. 

9. IF A is M and B is VH and C is H and D is H THEN E is VH. 

10. IF A is VL and B is L and C is L and D is VH THEN E is VH. 

1 1. IF A is VL and B is VL and C is VL and D is H THEN E is M. 

The plot of actual versus predicted sulphur for one of five classes yielded highest R 


value of 0.94. 



2.1.6 Optimization by GA 


GA works with a population of solutions, and during the iterations it continually 
updates the population according to the set objective function, throwing away poorer 
solutions and replacing them with improved solutions GA is, no doubt, computer 
time and memory intensive but that is not a limitation at all with small optimization 
problems discussed in this work. In fact GA is now accepted as one of the most 
powerful techniques of evolutionary search and optimization in process metallurgy. 


2.2 Plant data 

The plant data used in this work relates to 400 ton torpedo ladles (for details see 
2 . 10 , 2 . 11 ]^ which desulphurization was carried out by injecting calcium carbide 
powder through a submerged twin-whole lance. Nitrogen (2 m^/min) was used as a 
carrier gas. A typical plant data set is shown in Table 2.1. Data for 229 heats was 
collected, and the range of variation of each parameter was as follows: treatment time 
(473-1719 ± 20 seconds), hot metal weight (250-379 ±5 tons), initial sulphur content 
(90-440 ppm ± 5 ppm), CaD injection rate 45.0-71.0 ±1.5 kg/min), and final sulphur 
content (13-180ppm ± 2 ppm). 



2.3 Results and discussion 


A typical set of data for a sample calculation by different models is given below. 



Symbol 

units 

Value 

Treatment time 

tlotal 

sec 

1808 

Weight of hot metal 

W^ 

ton 

334 

Initial sulphur content of bath 

Si 

mass % 

0.033 

Injection rate of CaC 2 

Ir 

kg/min 

51 

Final sulphur content of bath 

Sf 

mass % 

0.004 


The models described in sec 2.3 are applied to predict the final sulphur content of hot 
metal. 

2.3.1 Application of models 

(a) Application of regression model to predict the required mass of calcium 
carbide for desulphurization 

The regression model developed from the stastical analysis (Eqn 2.1) may be used to 
predict the required mass of calcium carbide needed for desulphurization, Qcao, to 
achieve a desired final sulphur content of the bath. 

The following model, known as Ballistic model, is used 

'' S 

QcaD = Xi + Wm X In -i- X {X 2 +X 3 xSi + (X4 + X5xlR)xSf+X6xlR} (2.1) 



Where, Xi Xg: Constants, Wm: Hot metal weight, S,; Initial sulphur content, Sf: 

Final sulphur content, Ir: Injection flow rate, Qcao: Mass of calcium carbide needed. 
For the data set used in the present work, the values of the estimated parameters 
obtained by multiple linear regression are: Xi = 103.12, X2 = 4.7 x lO"^, X3 = 0.014, 
X4= 0.076, X5 = -0.081, X6= 0.0012. 

Thus, the predicted value of the required mass of calcium carbide for 
desulphurization, Qcao, can be calculated by using Eqn (2.1). 

^ S '' 

QcaD = Xi+WmXln X {X2 + X3 X Si + (X4 + X5 X Ir) X Sf + Xg X Ir} 

= 103.12 + 334 X 10^ X In x {4.7 x lO"^ + 0.014 x 0.033 

U.004j 


+ (0.076 -0.081 X 0.85) x 0.004 + 0.0012 x 0.85} 

= 1492.95 kg 

From the data the actual required mass of calcium carbide for desulphurization, Qcao, 
is Actual, QcaD = tiotal X Ir 



x51= 1567.8 kg 


The plot of actual versus predicted mass of calcium carbide is shown in Fig. 2.2. It is 
possible to back calculate the final sulphur content of hot metal (Sf) from Eqn (2.1) 
and from the plot of actual versus predicted sulphur, the value of = 0.87 was 
obtained. This value is marginally higher than the value of R^ = 0.82, as 

determined in the earlier regression model (in which no logarithmic relationship was 
assumed between Si and Sf). It demonstrates that if the process fundamentals are 
properly understood then it is easier to select the regression parameters for improving 
the accuracy of prediction. 



(b) Application of the model based on second order chemical reaction to 


predict the final sulphur content of the bath 



From the above relationship the value of final sulphur content of the bath, Sf, is 
calculated by iteration, and Sf = 0.0037; the actual sulphur content observed in this 
heat was 0.004%, and the two values (actual and predicted) agree well. 

Results for all the heats are plotted in Fig. 2.4 and it can be seen that the overall 
value is 0.80. In this kinetic model no separate treatment is incorporated to account 
for the specific contributions of transitory reaction (while the particle is traveling in 
metal) and the permanent contact reaction (when the particle has joined the top slag 
reactions). 



Predicted Sr, % Predicted Qcad (Kg) 



Actual Qcad (kg) 


Fig. 2.2: Actual versus predicted Qcao by regression model 



Fig. 2.3: Actual versus predicted Sfby regression model 





(c) Application of the model based on mass transfer (first order), including 
the transitory and permanent reactions 


The model based on second order chemical reaction (Eqn 2.4) may use to predict the 
final sulphur content of the bath 


■Ml d[%S] 


M. 


dt 


-Rts +Rp +Rb 


Where, Rp is the contribution of individual powder particles while they are rising 
independently in the liquid metal, Rb is the contribution of the powder particles which 
are captured in the rising carrier gas bubbles, and Rts is the contribution of particles 
when they join the top slag. 

For a small time step At = 0.1 sec the values of Rts, Rp, Rb can be calculated by. 


M 100 


OS ~ - 


Ml([%S], -[%s]) + M,(.,0)(%S), 

L“(M,(t=o)+lRt) 



Fig. 2.4: comparison of the actual versus predicted sulphur for the second order 

[ 2 . 1 ] 


chemical kinetics model 




The list of various symbols and their values used in the above equations is given in 
Table 2.4. By putting the values in the above equation, 


„ 8.70x 10“'' X 40x7000 

Rts = X 

32x100 


0.033 • 


V 


334x10^ (0.033 -0.033) + 5x0 
2000 X (5 + 0.85x1) 


= 0.0392 


R, 




M .100 


1-exp 


r \ 

6kppLtrp 
L'dpPp 


The list of various symbols and their values used in the above equations is given in 
Table 2.4. By putting the values in the above equation, 

_ (1 - 0.3) X 0.85x2000 
32x100 


1-exp 


^ 6x2.36xl0-'x7000 ^ 


2000x15x10-' xl740 


X [0.033] 


= 0.045 


f-RL^s 


M.lOO 


1 - exp 


^ 2.38K,p,trbgr^Pp^^ 


V 


273d, RL^f R 


xi%; 


The list of various symbols and their values used in the above equations is given in 
Table 2.4. By putting the values in the above equation. 


R. = 


0.3x0.85x2000 


32x100 


1-exp 


2.38 X 2.36 X IQ-ROOO x 7.13 x 0.033 x 1473 
273 X 0.02 X 0.85 X 2000 x 0.3 




xl 


X [0.033] 


)) 


= 0.0053 

Finally, the decrease in sulphur content of metal small time step. At, is 


A[%S] = (R^+Rp+Rb)— At = (0.0392 + 0.0405 + 0)x- — ^^-^xO.Ol 

Mg 


334x10' 


= 3.92x10'® 

For the total time the predicted value of the final sulphur content of the bath, Sf, is 


0.039 % 



The results obtained for all the heats are plotted in Fig. 2.5. The value of = 0.84 
(obtained from the plot of actual versus predicted sulphur) is higher than that of the 
second order chemical kinetics model (R^ = .80) but lower than the two regression 
models discussed earlier. This first order mass transfer model, even though more 
sophisticated in approach to account for transitory and permanent contact reactions, is 
based on several assumptions and has some parameters which can be tuned or 
optimized. The advantages of optimization of this model by GA will be discussed 
later. The results of some other models based on first order mass transfer are 
discussed elsewhere 

(d) Application of ANN models to predict the final sulphur content of the 
bath 


Results of the three typical ANN models are presented in Fig. (2.6-2. 8) for three layer 
(4-5-1 net), three layer (4-6-1 net), and four layer (4-6-3-1 net), respectively. Simple 
back propagation method was used for training by using a sigmoidal function, with 
momentum factor a = 0.5, training rate p = 0.5, and logistic T = 1.0. The training set 
contained 91 data and test set contained 79 data. 

The R^ = 0.86 is obtained for the case of simple net 4-5-1. The other two nets also 
gave the R^ values close to 0.85. 



Predicted Sf, % Predicted Sf, 



Fig. 2.5: Actual versus predicted Sf by mass transfer model 



Fig. 2.6: Actual versus predicted Sfby ANN 4-5-1 net 





Predicted Sf, % Predicted Sr, % 



Fig. 2.7: Actual versus predicted Sfby ANN 4-6-1 net 



Fig. 2.8: Actual versus predicted Sfby ANN 4-6-3 -1 net 






(e) Application of Fuzzy rule based expert system (FREES) model to predict 


the final sulphur content of the bath 


The plot of actual versus predicted sulphur (for the data set given in the Table 2.1) is 
shown in Fig. 2.9. 



Fig. 2.9: Actual versus predicted Sf by fuzzy rule based expert system (FREES) 


model 




(f) Application of genetic algorithm (GA) to predict the final sulphur content 


of the bath 


In an earlier work the model based on mass transfer was optimized by using GA. 
A simple genetic algorithm (SGA) was employed with Roulette Wheel Selection, 
crossover probability ranging from 0.7-0.85, and mutation was less than 0.01. The 
relative contributions of powder particles (see equation (2.4)) rising independently in 
the liquid metal (Rp), powder particles captured in the rising carrier gas bubbles (Rb), 
and the powder particles when they join the top slag (Rts) were investigated. It was 
shown that the contribution of particles captured in rising gas bubbles should be 
neglected. This significant result could not be obtained by the application of a 
regression technique even after optimizing the coefficients Rp, Rb, and Rts in 
equation (2.4) The result predicted by the combined GA-kinetic model has been 
found to be in agreement with the experimental observations that there is a critical 
limit to the ratio of mass of powder to the mass of gas injected per unit time at which 
the utilization of desulphurization agent is optimal with respect to desulphurization 
time and the cost of reagent injected. In the final GA-optimized model the kinetic 
term Rb was dropped altogether The GA optimization further revealed that the 
contribution of top slag (Rts) was nearly 1.5 times of the contribution of particles 
independently rising in the metal (Rp). The optimization of the exponent of Rp and Rts 
improved the predictions further and from the plot of actual versus predicted sulphur 
the R^ value obtained was 0.80. The benefit of optimizing the mass transfer model by 
using GA was two fold: better understanding of the process kinetics and an improved 
kinetic model based on process fundamentals. 



2.4 Comparison of modeling strategies 


The results of both conventional and AI techniques demonstrated in this work are 
summarized in Table 2.3, in terms of average error, mean square error and the square 
of the product-moment correlation coefficient (R^) for actual versus predicted final 
sulphur content. To start with, the experimental data was selected for a wide range of 
final sulphur content, 0.005- 0.018 %. It was, however, observed that in all the 
models, in the range of final sulphur lying between 0.010-0.018 mass %, the predicted 
out put was almost independent of the actual out put. For example, from the results of 
multiple linear regression analysis of the two above mentioned data sets, for final 
sulphur content < 0.010 and 0.01-0.018 mass % respectively, a considerable decrease 
in standard error of estimate was observed for the first set but for the second set, 
0.01-0.018 mass %, the value dropped to 0.0217. This suggested that the process 
data should be split in two parts for prediction purposes namely, 0.005- 0.012% and 
0.01-0.018%, respectively. In the present work the results are presented for the former 
case when final sulphur content is less than 0.012%. It is advisable to tune or 
optimize a model to experimental data in specific ranges if that helps to increase the 
predictability. The process characteristics can change in different regions. 

The plus point of a kinetic model based on process fundamentals is that as the process 
understanding improves it becomes possible to make further modifications in the 
design/technology. For example, having learnt firom GA optimization experiments 
that the contribution of powder particles rising in the bubbles was nearly zero, the 
lance nozzle design was improved to inject powder deeper into the melt and also 
increase residence time of the particle freely rising in the melt. This led to a 
significant improvement in process efficiency. 



If we compare only the numerical values of for conventional models, both the 
multiple regression models (R^ = 0.87) and second order chemical kinetics model ((R^ 
= 0.80) have performed well. The process knowledge gained with these models can 
however be misleading at times. The same is true for AlSDSf model (R^ = 0.88). A 
major limitation of using ANN models in a plant environment is that in certain 
regions of the multidimensional space, wherein enough training data is not available, 
the predictions can become totally unreliable. The fuzzy rule based expert system 
model appears to give the best result (R^ = 0.94) because it was possible to make 
specific rules based on actual plant data and then focus on a small class. For example, 
in the case of the data in the range of 0.01-0.018 mass %S, the base of the triangular 
membership function for final sulphur was increased as the final sulphur content 
increased. In spite of better prediction of final sulphur content, FREES model does 
not provide any additional information or insight into the process. The rule base in 
FREES also needs to be updated frequently in a plant environment and this is not a 
simple exercise. From the point of view of practical convenience, therefore, only 
when no other suitable model is available, one should try the FREES approach. 
Eesides a high R^ value, the other selection criteria of a good model is that both 
average error and mean square error should be as small as possible (given in Table 
2.3), and the spread of error should be uniform over the entire range of predicted 
variable. The residual error plots for different models are compared in Fig. 2.12. The 
spread of error is almost uniform in the case of ANN model followed by GA-kinetic 
model. The second order chemical kinetics model shows the worst spread followed by 
Sheshadri/Deo kinetics model, ballistic model and the second order chemical kinetics 


model. 






Second Order Chemical Kinetics Model ANN 4-5-1 Net 



FRBES model GA-kinetics model 


Fig. 2.12: Plot of residual error (i.e., actual-predicted) versus heat number for 


different models 






Table 2.1: Typical plant data set on desulphurization 


Total no. of observations: 1 1 

Treatment 

Hot Metal 

Initial Sulphur 

CaD injection 

Final Sulphur 

time (Sec) 

Weight (Tons) 

(mass %) 

rate (kg/min) 

(mass %) 

1404 

280 

0.026 

67 

0.003 

1421 

309 

0.024 

70 

0.003 

473 

303 

0.016 

62 

0.005 

1161 

362 

0.030 

73 

0.005 , 

1085 

336 

0.028 

63 

0.008 

559 

308 

0.021 

63 

0.009 

483 

304 

0.018 

58 

0.011 

501 

317 

0.020 

63 

0.011 

779 

295 

0.030 

68 

0.012 

796 

344 

0.033 

60 

0.015 

__ 

1 

271 

0.022 

68 

0.017 










Table 2.2: Fuzzification details of treatment time 


^(VL) 

1 

(675-X)/202 

0 

X<473 

473 < X < 675 

X = 675 

^(L) 

(X-625)/67.5 

625 < X < 692.5 


1 

X = 692.5 


(760-X)/67.5 

692.5 < X < 760 

|a(M) 

(X-725)/137.5 

725 <X< 862.5 


1 

X = 862.5 


(1000-X)/137.5 

862.5 <X< 1000 

^(H) 

(X-950)/225 

950<X<1175 


1 

X=1175 


(1400-X)/225 

1175 <X< 1400 

li(VH) 

0 

X=1300 


(X-1300)/419 

1300 <X< 1719 


1 

X>1719 




Table 2.3: Comparison of the performance of different models 


MODELS 

Average 

Error 

Mean 

Square 

Error 

R" 

Ballastic model (Linear Regression) 

0.0016 

1.21 X 10'^ 

0.87 

Second order chemical kinetics 

0.0022 

1.77 X 10'^ 

0.80 

Seshadri/B.Deo, mass transfer 

0.0013 

9.98 X 10’^ 

0.84 

ANN [ 4-5-1] 

0.0044 

1.26 X 10'^ 

0.88 

ANN [ 4-5-3-1] 

0.0072 

1.6 Ox lO’’^ 

0.86 

ANN [ 4-6-1] 

0.0012 

1.01 X 10'^ 

0.81 

ANN [ 4-6-3- 1] 

0.0040 

1.27 X 10'^ 

0.80 

FREES model’ 



0.94 

Optimization by genetic algorithm 

0.0014 

1.09 X 10'" 

0.80 

* best result obtained for one of the five data set 

s (class) in FRl 

BES 





Table 2.4: The list of various symbols and their values used in equation (2.4) 


Symbol 


Value 

Af 

Area of slag metal interface (m'^) 

40 

dp 

Average particle diameter (m) 

1.5 X 10'" 

f 

Fraction of particles entrapped in the bubbles 

0.3 

h 

Depth of injection (m) 

1.8 

Ir 

Rate of injection of the desulfurizing agent (kg/s) 

0.85 

kb 

Sulfur partition coefficient between bubble 

2.36 X 10'^ 

kf 

Sulfur transfer coefficient through the metal slag 

8.70 X 10-^ 

kp 

Sulfur transfer coefficient at particle melt boundary layer 

2.36 X 10'^ 

Lg'’ 

Sulfur partition coefficient between particles and melt 

respectively. 

2000 

Ls® 

Sulfur partition coefficient between particles entrapped in 

the bubble and melt. 

2000 

J-^s 

Sulfur partition coefficient between top slag and melt 

2000 

McaC2 

Molecular mass of CaCa (g/mol) 

64 

Ml 

Weight of metal (Ton) 

334 

Ms 

Molecular mass of sulfur (g/mol) 

32 

Mts(t=o) 

Mass of carry over slag (kg) 

5 

Po 

Average pressure in side the bubble (atm) 

1.0 

Q 

Gas flow rate (Nm'^/s) 

0.033 

Sf 

Final Sulphur content (%wt) 

0.004 

[%S]o, 

Si 

Initial sulfur concentration in the melt (%wt) 

0.033 





t 

Total time for desulphurization treatment (s). 

500-3400 

PL 

Specific weight of melt (kg/m^) 

7000 

PP 

Specific weight of particle (kg/m^) 

1740 



Chapter - 3 


Theoretical and practical aspects of nitrogen 
and oxygen absorption, and inclusion 
engineering 

3.1 Model of Oxygen and Nitrogen Absorption by Liquid Steel 
during Tapping 

3.1.1 Introduction: Liquid steel, at the end of heat in a EOF, contains dissolved 

oxygen (400-900ppm), hydrogen (2-6 ppm), nitrogen (40-60 ppm) and the 
temperature of metal varies in the range 1 873-1 973K. During tapping of the liquid 
steel from EOF into a ladle, the falling stream of liquid steel is surrounded by the 
gaseous atmosphere (Fig. 3.1), which is rich in both Na and O 2 . There is diffusion of 
Na and Oa into the stream of liquid steel from the atmosphere and hence additional 
nitrogen and oxygen are picked up by the steel. Also, the metal bath in the ladle 
becomes turbulent during tapping due to of kinetic energy of falling stream and 
energy input due to bottom stirring and this leads to entrapment of gas bubbles in the 
metal. The gas bubbles may dissolve (partially or fully) and therefore further add to 
the dissolved nitrogen and oxygen content of steel. 

The phenomenon of Na and Oa pickup was studied in the past by several investigators 
[3.1-3.6] mathematical model of Na and O2 pickup presented in this work has been 
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adapted ' ' and further improved to incorporate the effect of absorption by tapping 

stream as well as absorption through the top surface of metal bath. 

3.1.2 Description of the kinetic model of absorption 

During tapping, the N2 and O2 dissolution in liquid steel occurs simultaneously by: 

a) Diffusion through the surface of the falling stream 

b) Diffusion through top surface of the molten metal pool in the ladle 

c) Diffusion of oxygen and nitrogen from air bubbles, which are entrained in the 

molten pool. 

Following assumptions have been made in the kinetic model developed in this work: 

• The liquid steel in the ladle is uniformly mixed during tapping due to high 
amount of turbulence 

• Temperature distribution is uniform throughout the system and almost 
constant during tapping 

• Tilting angle of the EOF vessel changes such that metal flow rate remains 
steady (constant stream diameter) throughout the tapping operation 

• All the oxygen in the entrapped bubble is completely absorbed, but nitrogen is 
absorbed only partially. 



Fig. 3.1 : Schematic diagram of tapping of liquid steel into a ladle 






Different mechanisms of absorption of nitrogen and oxygen have been discussed 
below 

(a) Absorption of nitrogen and oxygen by diffusion through the surface of the 
falling stream 


The equilibrium constant for oxygen dissolution in liquid steel is 
!4(02)g = 2[0] Ko = 



= 10 

11370 

-t- 3.645 

T 

P02 ) 




(3.1) 


Where, ICo is the equilibrium constant, Po^ is the partial pressure of oxygen in the 

ambient atmosphere, [0]' is the concentration of oxygen at the air- stream interface 
within a small height Ah, and T is the temperature of metal. Since the chemical 
reaction of oxygen dissolution is extremely fast the oxygen absorption by the 
liquid steel stream surface is assumed to be controlled by mass transfer of oxygen in 

the metal phase. From Eqn (3.1), if we know the values of Po^ and T, the interfacial 
concentration of oxygen at the metal surface can be calculated. The increase in 


oxygen content, A[<9]„ , during the time step At can be calculated from 

A[6],s=f;ko^([0]'-[0]‘)At (3.2) 

^ AVh 

Where, [O]'' is the concentration of oxygen in a small volume, AVh of height Ah of the 
liquid stream. Ah is the area of the small volume exposed to the atmosphere, and ko is 
the coefficient of mass transfer of oxygen in the stream of liquid metal. 



The equilibrium constant for nitrogen dissolution in liquid steel is 


N2(g) = 2[N] 


Kn = 


'Mill' 

Pn2 J 



371.9 

T 


- 2.5 


(3.3) 


Where, Kn is the equilibrium constant for the dissolution of nitrogen by liquid steel, 
is the partial pressure of nitrogen in the ambient atmosphere, [N]' is the 

concentration of nitrogen at the air-stream interface within a small height Ah. The 
phenomenon of absorption of nitrogen by the surface of metal stream, similar to that 
for oxygen, is also assumed to be controlled by mass transfer of nitrogen in the metal 

phase. From Eqn (3.3), for known values of and T, the equilibrium concentration 

of nitrogen at the metal surface can be calculated and the rate of nitrogen absorption 
can therefore be calculated from 

A[N]Ts=i;k„^([N]‘-[N]'')At (3.4) 

Where, [N]^ is the concentration of nitrogen in a small volume, AVh of height Ah at 
the liquid stream. Ah is the area of the small volume exposed to the atmosphere, kw is 
the coefficient of mass transfer of nitrogen in the stream of liquid metal. 

The procedure adopted for oxygen and nitrogen pickup (flow chart given in Fig. 3.2) 
is as follows: 

Within a small time step At the increase in oxygen and nitrogen content of the stream 
of liquid metal for a small stream length. Ah, can be calculated from Eqns (3.2) and 
(3.4), respectively. The increase in [(9]’’ and [A]*’ within the small stream length. Ah, is 
continuously updated. 



(3.6) 


[Nf =[Nf 

The new values of [O^ and [A/]’’ , thus obtained, are used for calculating the [Of’ and 
[A^’’ values over the following small stream length ( i.e. Ah to 2Ah). This iteration 
process continues until the total stream length equals to the stream height, Zt, i.e. 
when the stream joins the metal pool of ladle. The total oxygen and nitrogen pickup 
due to absorption through the entire surface of the stream, over a small time step At is 
therefore, 


^^^Stream ~ 





(3.7) 




/i=0+A/ 


(3.8) 


The total stream height, Zt, changes as a function of time because of the filling up of 
the ladle and the change in tilt angle of the BOF vessel. The value of Zt, at any given 
time, is 

Zt = (LH-ZL) + ZF + Zc (3.9) 

Where, Lh is the height of the ladle, Zl is the height of the metal pool of the ladle, Zf 
is the distance between top of the ladle and the bottom of the converter, Zc is the 
distance between the position of the tap hole of the converter and the ladle (see Fig. 
3.1) at the horizontal state and that during tapping. The height of the metal pool of the 
ladle, Zl, can be calculated from 


^L|/+Ar 


AVr 
Z ,+ — ^ 

^ Ar 


t 


(3.10) 



Where, is the interfacial concentration of oxygen at the metal pool and the 

atmosphere, is the average bulk concentration of oxygen in the liquid bath, Al 

is the area of the surface of liquid pool, ko is the coefficient of mass transfer of 
oxygen at the pool of liquid steel. 

The equilibrium constant for oxygen dissolution is given by Eqn (3.3). The 
phenomenon of absorption of nitrogen by the surface of metal stream is also assumed 
to be controlled by mass transfer in the metal phase. From Eqn (3.3), for known 

values of and T, equilibrium concentration of nitrogen at the metal surface can be 
calculated and the rate of nitrogen absorption can be calculated jfrom 

A[N], = k, ^([N]i.„ - [N]“)At (3.12) 

Where, the concentration of nitrogen at the metal pool-air interface, [N]^^ is 

the average bulk concentration of nitrogen at the liquid bath, Al is the area of the 
surface of liquid pool, kn is the coefficient of mass transfer of nitrogen at the pool of 
liquid steel. 

(c) Absorption of nitrogen and oxygen from gas bubbles entrained in liquid steel 

The rate of gas entrainment by a stream of liquid steel into a pool of liquid steel, , 
as described by Choh et al is 

V° 

V = — ^ 

® sine' 


(3.13) 



Where, V° is the rate of the gas entrainment of the vertical component of the stream 

and 0^ is the incident angle of the tapping stream with respect to the liquid steel pool 
in the ladle. V ° can be calculated from 


V“ = 0.02V, 


R -a 


-f3 


(3.14) 


Where Rc radius of the cavity of the region in which the falling stream impinging on 
metal surface, and a are the radius of stream at the tap hole exit and at distance Z 
from the tap hole exit. can be calculated from 


R, = 


(l. 14 + 8.9x10“^ -1.067 


0.5 


(3.15) 


Where Ei is the kinetic energy of the stream of liquid just before it heats the liquid 
pool. 



(3.16) 


Where U is the velocity of the stream of liquid metal just before it heats the liquid 
pool. The mass of the protuberant portion of the stream, Wm, and can be calculated 
from. 


m = TtpX 


las — 



(3.17) 


Where, X is the wavelength, are the amplitude of the surface disturbance of 

stream at the distance Z and at the tap hole exit, respectively. The radius of stream at a 
distance Z from the tap hole exit can be calculated from 
a = (3.18) 

Where, ^ is the function of ratio of total energy to kinetic energy (^ > 1), Z and 
velocity of the stream at exit of the nozzle, Uo- Ihe value ot c can be calculated by 



The wave length of ripples on the surface of metal can be calculated from, 
i = 8.38 a (3.20) 

The velocity of the stream, U, of liquid metal just before it heats the liquid pool can 
be calculated from 


17 = 17 ,^^ 


(3.21) 


Amplitude of the surface disturbance of stream when the stream heats the metal pool, 
£ is given by 


£ = exp 




(^’ -l)xFrxWe 


0.35 




(3.22) 


Where, Fr is Froude number FV = , dn is the diameter of the nozzle at exit, g is 

the acceleration due to gravity. The weber number. We, is defined as 

We^HlPpA. 

Where, is the density of liquid steel, Cs-a is the surface tension of the liquid. 

Amplitude of the disturbance (ripple) on the surface of stream at the exit of nozzle, So, 
is given by 

(3.23) 

Where, Lq is Length of the nozzle, and Re is the Reynold’s number of the stream 


ao expl 


-44.3 


fr 

Vd„2 


Re 


- 0.1 


defined as 



Where, is viscosity of the liquid steel. 


The incident angle of tapping stream with respect to liquid steel pool in the ladle, 0^ , 
can be calculated from the following equation, 


0 = tan" 


(2gZ + U^ cos' 0/2 


Uq sin0 


(3.24) 


Where, as before, 0 is the exit angle of liquid metal stream through the nozzle. The 
increase in nitrogen content during the small time step At is. 


TlN,[%Nr]Pa.V, 0.0144r,^,V 


A[N],„ = 


W, 


- g.At= ^^At 


W, 


(3.25) 


Similarly, the increase in oxygen content during a small time step At can be calculated 
from 


a[6] 


T,‘J%or]p.tV, 0.00437Ti‘y, 


bbl 


■At = 




w, 


At 


(3.26) 


With the help of relationships (3.13)-(3.26) it is possible to calculate N 2 and O 2 pickup 
due to the bubbles entrained in the metal. The details of calculations are provided in 
section 3.5. 


(d) Total increase in dissolved nitrogen and oxygen content of the bath 

The total increase in nitrogen content during the small time step At can be calculated 
from 

A[N], = A[N]ts + A[N]s + A[N]bb, (3-27) 

Hence, the total increase in nitrogen content of the liquid metal during tapping is 



(3.28) 


A[N] = Xa[N], 

t=0 

Similarly, the total increase in oxygen content during the small time step At can be 
calculated from 

A[6], = A[0]ts + A[6]s + A[6],y (3.29) 

Hence, the total increase in oxygen content of the liquid metal during tapping is 

A[6] = Xa[0], (3.30) 

t=0 

The different symbols, used for calculation of oxygen and nitrogen pickup, are listed 
in Table 3.1. The symbols are in CGS units as discussed in the original work 


Table 3.1: List of symbols used for calculation in the absorption model 


Symbol 

Meaning 

ao. a 

Radius of stream at the tap hole exit and at the distance Z from the 
tap hole exit, cm 

Ah 

Area of the stream of liquid metal at a small length Ah at height h, 
cm^ 

dk 

Average diameter of the entrained bubbles, cm 

dft 

Diameter of the tap hole, cm 

F, 

Area of the Gas-metal interface (cm"*) 

Fr 

Froude number 

Gt 

Total gas volume entrained, cm^ 

h 

Distance from the surface of liquid steel in the converter to the inlet 

of tap hole, cm 

H 

Height of the converter, cm 

Kg 

Mass transfer coefficient in gas phase, cm/sec 

m-N 

Mass transfer coefficient of nitrogen in the liquid phase, cm/sec 

mo 

Mass transfer coefficient of oxygen in the liquid phase, cm/sec 

1 

Distance in the longitudinal direction from the bottom of the 

converter, cm 

lo 

Length of the tap hole, cm 

m 

Mass of the protuberant portion of the stream, g 

[%Nf] 

Nitrogen present in the air by mass percent, % 

m m 

No^ , N N 2 

Rates of mass transfer in the liquid phase and the chemical reaction 

at the interface, respectively, moiy(cm .sec) 

a[N],„a[nu 

Average nitrogen absorbed from the entrained gas bubbles and 
through the surface of the tapping stream, % 

[0]\[0]‘ 

Oxygen content in bulk of liquid steel in falling stream and at the 
metal-atmosphere interface, respectively, % 

[%Of] 

Oxygen present in the air by mass percent, % 

A[d]„,A[G],„ 

Average oxygen absorbed from the entrained gas bubbles and 
through the surface of the tapping stream, % 



Pqi > Ki 

Oxygen partial pressure and its value at the interface 

Rc 

Radius of the cavity, cm 

Re 

Reynolds number 

Rl 

Radius of the ladle, cm 

T 

Temperature of the liquid steel, K 

Uo, u 

Flow velocity at the tap hole exit and its value at the distance Z 
from the tap hole exit, cm/sec 

V 

^ g ^ ^ g 

Rates of the gas entrainment of the vertical stream and of the tilting 
stream, respectively, cmVsec 

Vh 

Volume of the stream of liquid metal at a small length Ah, at a 

height h 

Vi 

Rate of tapping, cm^/sec 

We 

Weber number 

Wt 

Initial mass of the liquid steel in the converter, g 

Wl 

Mass of the liquid steel in the ladle at time t, g 


Weight of the liquid steel in the converter at time t, g 

Zt, Zf 

Tapping height and it’s final value, cm 

0,0^ 

the angle of BOF with the horizontal axis and incident angle of the 

tapping stream to the liquid steel pool in the ladle, degrees 

b b 

Absorption efficiency at the bubble-liquid interface of oxygen and 

nitrogen respectively 


Amplitude of disturbance present on surface of stream at a distance 

Z and at the tap hole exit, respectively, cm 

A 

Wave length of the ripples at the surface of liquid steel, cm 

^s-a 

Surface tension at air and liquid steel interface, dyn/cm 

pFe 

Density of liquid steel, g/cm 




3.2 Thermodynamics of deoxidation by aluminum 


The amount of aluminum reacted to obtain a given dissolved oxygen content, at a 
particular temperature can be calculated as follows 



2[Al] 

+ 3[0] 

-> [AbOa] 

(3.31) 




Let the initial cone: 

X 

Y 

0 

Let the final cone: 

X-2x 

r* • 

Y-3x 

X 


Let the initial moles of iron present be Z. The total number of moles of iron, 
aluminum and oxygen at equilibrium is Z+X+Y-5x 

The equilibrium concentrations of A1 and O in bulk are related to the equilibrium 
constants, Ke, 


K. =• 


‘■AUO, 


(h^)^(ho)^ 


(3.32) 


Where logK^ 


58473 

T 


17.74, 


(3.33) 


and h^i = f^, (mass% Al), ho = fo(inass%0) 

The equilibrium mass percentage of aluminum, [Al]e, can be calculated from 


[Al]e = 


(A -2x)* 27 *100 

27*(A-2x) + 16*(F-3x) + 56*Z 


(3.34) 


The equilibrium mass percentage of oxygen, [0]e, can be calculated from 


[0]e = 


(7 -3x)* 16 *100 

27*(A-2x) + 16*(7-3x) + 56*Z 


(3.35) 



On substitution to Eqn (3.32), 


f^*27*(X-2x)*100 


27*(X-2x)+16*(Y-3x)+56*Zj I 27*(X-2x)+16*(Y-3x)+56*Z 


fo*16*(Y-3x)*100 


(3.36) 


In our calculations activity coefficients of aluminum, (f^,), and oxygen, (fo ), and the 

activity of alumina, (a-AhOi ), are assumed to be unity for the shake of simplicity. The 
non-linear Equation (3.36) can be solved for one unknown x since the values, K^, X, 
Y, Z, fAi, fo are already known. 

The decrease in oxygen content of the bath due to deoxidation, AOt, is 
W 

A0T = -([0]i-[0]J^kg (3.37) 

Where, [0]i and [0]e are the initial and equilibrium mass percent of oxygen in the 
bath, and Wm is the weight of the liquid steel. 

The increase in aluminum content in the bath is AAlx 


W 

AAlT=([Al]i-[Al]J^ kg (3.38) 

Where, [Al]i and [Al]e are the initial and equilibrium mass percent of aluminum in the 
bath. 

Now, theoretically 1.125 kg aluminum is needed for removal of 1kg oxygen. Hence, 
theoretical aluminum needed for deoxidation, Tai, is 

Tai =1 .125 |AOt| + AAIt , (3-39) 

In a practical situation, a part of aluminum charged into the steel may bums or react 
with the atmosphere. If the efficiency of aluminum utilization is Pal then the amount 
of aluminum needed to be charged is 

A1 = n T (3-40) 

spT#5Tr«r 


Aluminum, the deoxidizer, is usually charged into the liquid steel during tapping in 
the form of aluminum ingots (20 kg blocks, at Bokaro steel plant). The aluminum 
blocks dissolve rapidly into the steel due to the turbulence created by the tapping 
stream. The calculations of efficiency of aluminum utilization (recovery) for the 
situation of Bokaro Steel Plant are discussed in sec 3.5.3 

3.3 Nucleation and growth of alumina particles in the melt 

The deoxidation reaction (of dissolved aluminum with dissolved oxygen) proceeds 
through the following steps 

1 . Nucleation of the solid almnina particles 

2. Subsequent growth of inclusions by mass transfer 

3. Growth of inclusions by collision with other inclusions (Brownian motion, Stokes’ 
collision, laminar shear etc.) 

4. Assimilation of inclusions by top slag 

Nucleation of alumina particles is discussed first (section 3.3.1), followed by 
development of the kinetic model and its application to 300 ton ladle (section 3.5). 
The phenomenon of growth and also removal of inclusions is reviewed in section 
3. 3. 2-3 .4.4, followed by application of different models in section 3.5. 


3.3.1 Nucleation of solid alumina particles in molten steel 


According to the classical theory of homogeneous nucleation, the solid alumina 
inclusions may start to nucleate from the parent phase (melt) when the change in the 
Gibbs free energy for the deoxidation reaction (Eqn 3.31) is negative. At equilibrium, 


dAG 

dr 


0 


(3.41) 


Where, r is the radius of the nucleous. The shape of the nucleating particles may be 
spherical, planar, cylindrical or star shaped. For the sake of simplicity, we assume the 
shape of the nuclei to be spherical. In the case of spherical particles, the sum of 
volume of free energy and the surface energy created owing to the formation of new 
surfaces is given by, 


AG =—7Tr^AG^+Aw'cr (3.42) 

3 

Where, AGv is the change in free energy per unit volume, and Os-i is the surface 
tension between AI 2 O 3 particle and liquid steel. 

On solving equation (3.41) and (3.42) 

AG,=-^ (3.43) 

Tc 

The value of AGv is given by; 


AG=-— Inn (3-44) 

y. 

The super saturation, H, is the ratio of the actual and the equilibrium concentrations 



Assuming that AI 2 O 3 is pure, aAuos = 1 • Thus, 


C. K _ ([mass%Al],)-([mass%0],)^ 

Ce, K, K, 

Where, [mass%Al]j, is the mass percentage of aluminum in bulk, [mass%0]j,is the 

mass percentage of oxygen in bulk. The value of Ke can be calculated from Eqn 
(3.33). From equations (3.44), (3.45) and (3.46), the following expression is obtained 
for the critical radius, rc ,of the nucleous: 


RTlnn 


RTln 


([mass% Al]b )^ ([mass%0]b )^ 


(3.47) 


The critical radius is defined as the minimum radius of a nucleus for which the 


thermodynamic conditions are satisfied, and, therefore, nucleation can occur only 
when r > rc. From equation (3.47) we can easily see that the critical radius of a 
nucleus decreases with increase of super saturation and with decrease of surface 
tension. Also, since super saturation changes both with time and temperature during 
the deoxidation process, the critical radius also changes accordingly. 

The value of the critical radius is the starting radius of each inclusion. Subsequent 
increment to radius (i.e. growth of nucleous) takes place by mass transfer. 


3.3.2 Growth of nucleous: 

Once a stable nucleous has formed, it begins to grow by diffusion of [Al] and [O] to 

AI2O3 

2[Al]+3[0] ^ [AI 2 O 3 ] 



Aluminum and oxygen dissolved in the melt combine at the interface of the inclusion 
to form AI 2 O 3 . A boundary layer develops at the surface of the inclusion and metal, 
as shown in the diagram (Fig. 3.2). Diffusion of any one or both the elements (A1 and 
O) across the boundary layer may be considered as the rate determining step. 


Inclusion 
Boundary layer 



Inclusion/layer interface 

Bulk metal/ layer 
interface 


Fig. 3.2: The boundary layer formed at the surface of the inclusion and metal 


3.3.3 Calculation of interfacial concentrations of oxygen and aluminum 

(a) Calculation of interfacial concentration oxygen 

Similar to the method discussed at section 3.2 the equilibrium mass percentages of 
aluminum and oxygen can be calculated from Eqns (3.34) and (3.35). 

The equilibrium constant, IQ, of reaction (3.31) from the Eqn (3.33). At the inclusion- 
liquid steel interface the value of IC is given by 



(3.52) 


Jo=ko(C^-Q) 


Flux balance of the mixed mass transfer controlled reaction is 

1 d[Al] ^ 1 d[0]_ d[Ah03] 

2 dt 3 dt dt 


(3-53) 


From Eqns (3.51-3.53), the equation for the interfacial concentration of aluminum at 
the particle-bulk interface, C^, , can obtained 
3.k^,.C^^,-2ko(C^-C^) 


r' = 

'^Al 


3.k 


(3.54) 


A! 


Where, C^, is the concentration of aluminum in bulk, Cq is the concentration of 
oxygen at interface, Cq is the concentration of oxygen in bulk, kAi, ko are mass 
transfer coefficients of aluminum and oxygen, respectively. 

From Eqn (3.54) the value of C^, can be calculated for a particular time step if the 
values of C^i , Cq ? Cq , kAi, ko are known. The value of C’^, and Cq can be calculated 

by iterative method by using Eqns (3.50) and (3.54). 

A fimction formed by using Eqn (3.54) is 

r:,r> 3.k„C>,-2k,(C>|-Ci,) 
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(3.55) 


The function F(Ca,) (Eqn 3.55) can be solved by using different starting values of 
C',, At the correct value of interfacial concentration the value of the function would 

Al ’ 



be zero. The corresponding value of C'q can be calculated by using C'^, , obtained 
from the iteration, in Eqn (3.50). 

The bulk composition of dissolved aluminum is given by 

Cl, = [mass%Al‘’ ]-^ (3.56) 

lOOxM^, ^ ^ 

Where, Mal is the molecular weight of aluminum. 

The composition of dissolved oxygen in the bulk is given by 

Cl = [mass%0^] ^ - ^ - gg; (3.57) 

100 xMq 

It can be seen from Eqn (3.54) that the interfacial concentrations of aluminxma and 
oxygen will continuously change with time as the bulk concentrations change 
continuously. 

3.3.4 Increase in radius of inclusion 


The total initial volume of an inclusion of radius r is 

Vol = — 7i*r^ *n 
3 

Where, n = number of inclusions present in the liquid steel. 

The increase in radius of an inclusion can be calculated as follows. 


(3.58) 
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Density of inclusion, p^, q , is 4000 kg/m^. Therefore, 


dr 51 


dt 4000 


kAi(C^A,-C’,) 


dt 


Or, ^ = 0.01275k^.(C^^,-C^) 


(3.59) 


From equation (3.59), the increment in radius, dr, can be calculated during each time 
step dt. After each time step, the bulk concentrations of oxygen and aluminum are 
updated. Growth of nuclei will occur untill equilibrium concentration is achieved in 
bulk. 


3.3.5 Change in bulk concentration of aluminum and oxygen in the bath 


From the flux balance the change in aluminum and oxygen in a small time step dt can 
be calculated by 

d[Al,03] _ 1 d[Al] d[Al303] _ 1 d[0] 
dt ” 2 dt ^ dt 3 dt 

Change in mass of alumina in the bath in a small time step dt is 

d[Al 303 ] = No^7t{[r + drf -r^}PA,203 (3-60) 



Where, No is the number of inclusion present at time t. r and (r + dr) are the radius of 
particles at time t and (t + dt), respectively. 


d[mass%Al] = ^ x 1 00 (3 6i) 

d[mass%0] = xlOO (3.62) 

Wm 

The change in bulk concentrations of aluminum and oxygen is updated and can be use 
for the next time step. 


3.4 Growth of inclusion by collision mechanisms 

As mentioned above, the size of nuclei increases initially due to the diffusion of 
aluminum and oxygen towards the nucleus-steel interface (as discussed in Section 
3.3.4). The inclusions can also grow in size by simultaneous contribution of one or 
more of the following mechanisms, depending upon the flow/mixing/stirring 
conditions present in the ladle ^ J 

Ostwald ripening (diffusion coalescence) 

- Brownian motion 

- Stokes collision 

- Turbulent collision 


Laminar shear 



3.4.1 Ostvvald ripening 


In Ostwaki ripening, small inclusions, which form in a system initially, slowly 
disappear because they combine (by diffusion coalescence) with larger crystals. The 
smaller crystals thus act as "nutrients" for the bigger crystals. As the larger crystals 
grow, the area around them is depleted of smaller ones. It has been reported that the 
relative contribution of Ostwald ripening to inclusion growth is almost negligible 
because it is a very slow process (atomic diffusion over short distances). 


3.4.2 Growth of inclusions by collisions caused by Brownian motion 


Collision of inclusions due to Brownian motion has been described by the Stoke- 
Einstein equation and the collision volume is given by 



(3.63) 


Where, k is the Boltzmann constant, and T is the temperature. In principle, the 
contribution of Brownian motion to inclusion growth should be limited to small 


particles (of approximately 1 micron) because only they can move randomly in a 
fluid. Further, for the contribution of Brownian motion to inclusion grovAh to be 
significant, the difference in inclusion sizes must also be large. For example, if p = rj 
in Eqn (3.63) the value of Wy increases only by a factor of 4 but if p = lOOrj then the 


value of Wy increases by factor of 100. Collision volume for Brownian motion nearly 
attains a constant value after approximately 1 pm (Fig. 3.3) 



3.4.3 Growth of inclusions by Stokes’ collision 


The inclusions rise through the steel bath because of the difference in density between 
steel and the oxide (alumina). The rising velocity (governed by Stokes’ law) increases 
with increasing size. Large inclusions rise at a faster rate than the small inclusions. 
Thus, during their rise large inclusions can collide with the smaller inclusions (called 
Stokes’ collision) and grow in size. 



Fig. 3.3: Collision volumes of different mechanisms (for an inclusion with a 2p.m 
radius) as a function of radius of the other inclusions 




The collision volume of inclusions formed by Stokes’ collisions, as suggested by 
Lindborg and Torssell 
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(3.64) 


Where g is the gravitational acceleration, g is the dynamic viscosity, and pfe and pi are 
the densities of the steel and the inclusion, respectively. It has been reported that for 
an inclusion of 2pm radius the contribution of Stokes’ collision to inclusion growth is 
significant only for inclusions of size grater than 5pm. The variation of collision 
volume with inclusion size is shown in Fig. 3.3; it is clear that only after 5pm microns 
the collision volume increases 

3.4.4 Growth of inclusions by turbulent collision 


The inclusions can collide with each other due to the influence of the turbulent eddies 
present in liquid bath. Saffman and Turner described a model for growth of rain 
drops by turbulent collision. The same model has been adopted for the case of 
inclusions in liquid steel The equation for the collision volume created due to 

turbulent collision is 


Wy = 1 .37r^a{ri O’ - 


^Fe 


(3.65) 



Where s is the turbulent energy dissipation, Vpe is the kinematic viscosity of steel, and 

a is the collision efficiency. In literature, a is reported to var\' between 0.27-0.65 for 
alumina inclusions in liquid steel. The effect of turbulent collisions is found to be 
most prominent on inclusions which are above 1-2 pm in size. 

3.4.5 Growth of inclusions by laminar shear 


The inclusions can move along with the different laminar streams in the metal bath. 
The laminar streams may form by stirring, due to temperature gradient or due to metal 
flow during teeming. An inclusion moving in a fast streamline can collide with an 
inclusion moving in a slower streamline. The collision volume of inclusions formed 
by collision due to laminar shear, Wy , as suggested by Smoluchowski is given by 



(3.66) 


Where, — is the velocity gradient in the direction where laminar flow exists. 
dx 


3.4.6 Decrease in the number of inclusions due to collision 

Let the initial inclusion size, after the attainment of deoxidation equilibrium, be ro and 
the number of inclusions present in the metal be n^,. Now, some of the inclusions can 
collide with each other (collision mechanism) and form bigger sized particles. If two 
no sized particles, in a small time step At, collide and form a new inclusion of size n 

then, r, = +rl = V2.r, 


(3.67) 



Again, after a small time step At, some of the newly formed particles (of radius ri) can 
collide with other inclusions present in the bath. Tf particles of radius ri and ro collide 
with each other then the resultant radius would be Vs.r^ . For the collision between 

two particles of same diameter, ri the resultant radius would beV4.r„, and so on. The 
general expression for rn can, therefore, be written as 

r„=Vm.r, (3.68) 

where, n and m are integers. The growth process may continue in this manner till 
solidification. The total number of inclusions present in the bath decreases due to 
merging of inclusions by collision mechanisms, and the rate of decrease of the 
number of particles can be calculated by the general equation 


dn 


(3.69) 


Where N, and Nj are the number of inclusions, respectively, of two different sizes, and 
Wjj^ is the collision volume of the particular mechanism (such as stokes, turbulent 
etc.). 


3.4.7 Decrease in the number of inclusions due to assimilation by top slag 


When an inclusion reaches the top surface of metal it can get assimilated into top slag. 
For a given bath height H, the decrease in the number of inclusions due to 
assimilation can be calculated from 





(3.70) 


It may be noted that the processes of decrease in number of inclusions by collision 
and due to assimilation by top slag go on concurrently. 



3.4.8 Growth of inclusions by simultaneous action of two or more mechanisms 


Different mechanisms of inclusion growth, as described above, can occur 
simultaneously. Lindborg and Torssell Model suggested in 1968 a model in which 
only Stokes and gradient collision were considered to act simultaneously. Iyengar and 
Philbrook^^ ®^ suggested, in 1972, a model in which Stokes and turbulent collision 
were considered simultaneously. In fact, since 1970, there have been many studies to 
predict the growth and removal of inclusions in liquid steel. These studies^^ ’^'^ ’®^ are 
not limited to the case of a ladle stirred by natural convection, but also extended to 
removal of inclusion to the slag, removal of inclusion to the refractory surfaces, 
removal of inclusion by gas purging (attachment of inclusions to bubbles: Bubble 
flotation models), removal of inclusion due to induction stirring, floatation of 
inclusions in tundish, etc. Most of these studies make use of the work by Lindborg - 
Torssell, and by Iyengar - Philbrook these two models are described in this section. 

(a) Lindborg and Torssell Model 

According to this model, the rate of change of radius of inclusion, before the 
attainment of deoxidation equilibrium, is obtained from flux equation based on mass 
transport control, 

dr 
dt 


D CJt)-Co 
r Cp-Co 


Where, 



Cx is the time-dependent solute concentration at a large distance from the particle, Co 
is the solute concentration adjacent to the particle, Cp is the concentration of solute 
inside the particle itself. 


Cco(0 = C„(0)-iV 


A.n 


r\t){Cp - Co) 


(3.72) 


The radius of the final particle, at equilibrium, can be calculated from 

, 1 / 




3 C^(0)-C, 


0 


A.TC.N Cp - Cq 


(3.73) 


Only Stokes and gradient collisions are considered in this model. The effect of 
turbulent eddies in the collision mechanism has been neglected. Also, the effect of 
Brownian motion has been neglected for the large size of inclusions. 

If the inclusions in the steel bath have a density function f(r) then the volume of the 
inclusions is 


1= ^f{r)^7c.r^dr (3-74) 

Now, according to Eqn (3.67), if two particles with radius ri and r 2 collide, the 
resultant size is R = , Because of collision there is a decrease in the density 

function f(ri) and f(r 2 ) and an increase in f(R). The rate of change in f(ri) will be the 
sum over all particles raidus rj colliding during time step dt 


df{rO 

dt 


-fin)\f{r2)W{r„r2)dr^ 

0 


(3.75) 


Where W(ri,r 2 ) is the collision volume which can be calculated for an assumed 
mechanism (Eqns 3.63-3.66). There is also formation of particle of radius ri due to 
collision of particles having radius smaller than ri. A single generalized formula has 
been obtained to calculate the change in density function as follows 



clfir) 

di 


7 r V °° 

1 ^ rfr, -/(r,) (r, ,r^dr, 

i> V 2 y n 


(3.76) 


Finally, a decrease in the value of f(r) takes place due to the assimilation of inclusions 
in the slag phase. The rate of change of f(r) due to this process is 
df{r) 


a, 


(3.77) 


Where, Ur is the upward velocity of inclusions according to Stokes’ law and h is the 
height of the melt. The overall value of rate of change of f(r) with time is obtained by 
adding Eqns (3.76) and (3.77). 


(b) Iyengar and Philbrook Model^^ *7 


The motion of inclusions in a laminar flow due to natural convection in the r and z 
directions of a cylindrical coordinate system is given by 

m-^^ = 6.7t.p.r.(u - V ) 
dt ^ ^ r rJ 

ni^ = (m-m0.g + 6.7T.p.r.(u, -vj 
dt 

Here z axis is in the vertical direction and r is radial direction from the centre of the 
ladle. 

If Uz and Ur are the velocity of the inclusion int the upward and radial directions, 
respectively, then the displacement of an inclusion particle from its initial position can 
be calculated from the following equations, 



dr , dz 

— = u,. and — = u, 
dt 7 dt 



According to Snioluchowski^' if Brownian motion and fluid dynamic interaction 
between the particles are neglected then the probability of collision between the r, and 
rm sized particle is 

p(i,m) = -jtT.G.(ri +r„)^N(r„) =VgrN(rm) (3.78) 

The particle collision frequency, p(i,m), due to turbulent flow field inside the steel 
bathf^ '°^ is 

4 f 

p(i,m) = -TC. - (r^ +r^)'N(r„) (3.79) 

3 \vj 

If there are N(ri) inclusions per unit volume with radius r,- and N(rm) inclusions per 
unit volume with radius r™, then according to Smoluchoski^^'®^ , the number of 
collisions per second between these two species of inclusions is, 

C(i,m) = VgrN(rn,)N(ri) 

If the collision efficiency in a convective flow is Egr (value less than unity), then we 
can modify the above equation as 

C(i,m) = Vgr N(rn,) N(ri) Egr (3.80) 

Where Vgr is the collision volume per unit time for an inclusion with radius q and rm 
under a velocity gradient 
The contribution to Vgr comes from two sources: 

1 . Collision in the peripheral or annular region (Gp) 

2. Collision in the axial or pipe region (Ga) 

Therefore, V^^ =-j7i:(ri +r,„)^(Gp +GJ (3.81) 

And C(i,m) = N(ri)N(r,,)|rt(r. +r,,)^(Gp +G JE^ 

Adding the contribution of Stokes’ collision to the above. 



The values of Egr ^ and Estk can be calculated from the following equations. 


E.= 


0.75xln(2.Stk)-'^ T ; 
1 + ^ + 1 + - 


V 


Stk- 1.214 




j 


V 


J 


1 + - 


+- 


1 


h ) 


" r 
1 + ^ 


2 


V 


r. 


(3.83) 


1 J 


Where rm < ti and. 



(3.84) 


Finally, the expression for rate of change in the number of inclusions with radius 


dN(ri) _ 1 

dt “2^ 

% 


fj Tj 

J jC(l,mi)dr^dr, 


- Jc(i,mi)dr^-f(s) 

0 


(3.85) 


The meanings of different symbols, used in this section, are described in Table. 3.2. 



Table 3.2: List of Symbols used in section 3.5.8 


Symbol 

Meaning 


Symbol 

Meaning 

C(i,m) 

Number of collisions per unit 

volume per unit time between 

inclusions with radius p and r^ 


Tmax 

Radius of the largest inclusion in 

steel 

Egr 

Collision efficiency due to a 

velocity gradient in liquid steel 


Re 

Reynolds number for the larger 

inclusion 

Estk 

Collision efficiency due to 

Stokesian rise 


Stk 

r., , , PfRe 

Stokes number = 

9.p 

G 

Average velocity gradient 


t 

Time 

Ga, Gp 

Average velocity gradient in 

axial and peripheral region 



Velocity of an inclusion in r and z 

direction 

g 

Acceleration due to gravity 


■ 

Average fluid velocities in r and z 

directions 

ks 

Constant in Stokes’ equation 

! 

i 

V 

Volume of an inclusion 

L 

Height of ladle 


Vgr 

Collision volume per unit time for 

an inclusion with radius rj and rm 

imder a velocity gradient 

mi 

Mass of inclusion 


Vcr 

Volume of the ladle 


Mass of liquid steel displaced by 

an inclusion 


r 

Radial direction 

N(r^) 

Number of inclusions with 

radius rm per unit volume of 

liquid steel 


z 

Vertical direction 






p(i,m) 

Frequency of collisions per unit 

time of an inclusion with radius 

Tj with inclusion r,,, 


£e 

Energy dissipated per unit time 

per unit mass 

Ap/L 

! 

Pressure gradient 


P, Ps 

Density of liquid steel and 

inclusion, respectively 

r 

Radius of an inclusion r 


P 

Viscosity of steel 

ti 

Radius of inclusion i; i = i, m, 1 


V 

Kinematic viscosity of steel 



3.5 Results and Discussion; Application to 300 ton BOF at 
Bokaro Steel Plant 

3.5.1 Prediction of nitrogen and oxygen pickup from ambient atmosphere during 
the tapping operation 


A user friendly computer program has been developed to calculate the pickup of 
nitrogen and oxygen by molten metal stream during tapping. Flow chart of the 
algorithm is shown in Fig. 3.4 and various steps of calculation are as follows 

Step 1: Enter relevant input through user fiiendly interface, and save to input file 
“Data.txt 

Step 2: Set initial time, to, is zero i.e., t = 0 and the time step At is 0.01 sec (say). 

Step 3: Calculate the total bath hei^t Zt at time t by Eqn (3.9) 


Step 4: Set the distance from tap hole, h, of a volume, AV of liquid metal 
corresponding to small height, Ah, as zero, i.e., h =0 

Step 5: Calculate the interfacial nitrogen and oxygen concentrations i.e. 
concentrations at the interface of air-tapping stream, by using Eqns (3.1) and (3.3) 

Step 6: Calculate the increase in dissolved nitrogen and oxygen content of the small 

volume (A[N]ts and A[6]xs, respectively), due to absorption through the surface of 
the tapping stream, in time At by using Eqn (3.2) and (3.4) 


Step 7: Increase in the nitrogen and oxygen content of the bulk in AV in At is 
calculated by using Eqns (3.5) and (3.6) 


Step 8: The value of h is incremented by Ah. If h < Zt then go to Step 3, else continue 



Step 9: Calculate the total increase of nitrogen and oxygen ot steel by the tapping 
stream, from Eqns (3.7) and (3.8) 

Step 10: Calculate the nitrogen and oxygen concentrations at the interface of air- 
metal bath by using Eqns (3.1) and (3.3) 

Step 11: Calculate the increase in nitrogen and oxygen content of the metal bath 

( ^[N]s A[6 ]s , respectively) in time step At, due to absorption from atmosphere, 

by using Eqns (3.11) and (3.12) 

Step 12: Calculate the increase in nitrogen and oxygen content of the bath, due to 

dissolution of the entrained gas bubbles (A[N]by and A[6]j,y, respectively) in time 

step At, by using Eqns (3.25) and (3.26). The parameters used in Eqn (3.25) and 
(3.26) can be calculated from Eqns (3.13-3.24) 

Step 13: Calculate the total increase in nitrogen and oxygen content of the bath during 
the small time step At by using Eqns (3.27-3.29) 

Step 14: The value of t is now incremented by At. If t < total tapping time then go to 
Step 3, else continue 

Step 15: Calculate increase in total nitrogen and oxygen content of the bath (A[N] 
and A[0] , respectively) during tapping period by using Eqns (3.28) and (3.30) 

Step 16: Display the results of total nitrogen and oxygen pickup within the tapping 
time in the output window 





Fig. 3.4: Flow chart of oxygen and nitrogen absorption during tapping of liquid steel 















It is assumed that all the nitrogen and oxygen gas present in the bubble is dissolved in 
the bath (i.e. 100 /o dissolution of the gas into the bath) and the other 
parameters/conditions (assumed in case of 300 ton BOF) are given below 


Symbol 

Meaning of the symbol 

Numerical 

value 

dn 

Diameter of the exit nozzle of the converter 

14 cm 

h 

Height of steel above the inlet of tap hole in the converter 

100 cm 

H 

Height of the converter 

600 cm 

Kn 

Mass transfer coefficient of nitrogen in the liquid phase 

0.0025 cm/s 

Ko 

Mass transfer coefficient of oxygen in the liquid phase 

0.0025 cm/s 

lo 

Length of the exit nozzle of the converter 

100 cm 

T 

Temperature of liquid steel 

1923 K 

tj 

Total tapping time 

12 min 

Zf 

The distance between top of the ladle and bottom of 

converter 

400 cm 

Zl 

The height of metal pool in the ladle 

600 cm 

Zt 

tapping height 

1479 cm 

^initial 

Initial Angle of ejection of the liquid metal throu^ the 

nozzle exit of the converter. 0 decreases linearly with 

tapping time 

15° 

7 

Pfb 

Density of liquid steel 

7 gm/cm 

Pfb 

Viscosity of liquid steel 

0.05 g/cm.s 

^s-a 

Surface tension between air and liquid steel interface 

5 dyn/cm 


At 


Small time step 


0.1 sec 


(a) Magnitude of nitrogen and oxygen pickup through the surface of falling 
stream 


The increase in oxygen content, A[0]2.j , during the time step At can be calculated 
from Eqn (3.2) 

A[0]Ts = i;Ko^{[0]‘-[0]‘')At 

Where, Ah and Vh are the exposed surface area and volume, respectively, of a small 
length of the tapping stream. Ah. The value of Ah is calculated by dividing the total 
stream length into 500 equal segments (decided a priori) 



The value of Zt for the first time step. At = 0.1 sec, can be calculated from Eqn (3.9) 

Zt = (Lh - Zl) + Zf + Zc 

Where, Lh is the height of the ladle, Zl is the height of the metal pool of the ladle, Zf 
is the distance between top of the ladle and the bottom of the converter, Zc s the 
distance between the position of the tap hole of the converter and the ladle 
Now, Lh = 600, Zl= 0, Zf = 400, Zc = H cos 0 -Lo, H=600, 9 = 15°,Lo =100 
Hence, Zt = (600 - 0) +400 +[600 cos (15°) -100] 

= 1479.56 cm 

Therefore, the value of Ah is 2.96 cm 

The interfacial mass % of oxygen can be calculated from Eqn (3.1) 


["([onAjo 

[j'”0 3.645] 

T 

V Po2 / 




r , 3.645I 

OraOj^ -- VPo- • 1" ' ' 

Since the temperature of steel T = 1923K and po 2 = 0.2 1 , hence 



= 0.0337 % 

The initial mass % of oxygen in the bulk is 0.07% 

The increase in oxygen content, A[d]j^ , during the time step At for a small length of 
the tapping stream, Ah, is 

A[6],3=ko^([0]'-[01-)At 

= ko-^([0]'-[0]'')At 

Ah 

Since [O]’ = 0.0337, [O]’’=0.07, we get 
A[0 ]ts =0.0025^(0.0337-0.07)x0.1 
= 3.14xl0'^/„ 

Similarly, the increase in the nitrogen content, A[N]ts , during the time step At can 
be calculated from Eqn (3.4) 

A[N],3=t;k,^([Nr-[N]‘)At 

AV, 

The interfacial mass % of oxygen can be calculated from Eqn (3.3) 




Since in air pN 2 = 0.79, we get 


[N]' = 

= 0.0405 % 

The initial mass % of nitrogen in the bulk is 0.005%. Hence the increase in nitrogen 

content, A[N].j. 5 , during the time step At for a small length of the tapping stream, Ah, 
is 

A[N],s=k„^([N]‘-[N]‘)At 

= k„^([N]'-[Nf)At 
Ah 

On substitution of the values calculated above, 

A[N]ts = 0.0025 ^(0.0405- 0.005) X 0.1 
= 2.99x10-®% 

A typical profile of the increase in oxygen content caused by absorption through the 
surface of the tapping stream versus the distance from the exit of the tap hole of 
converter is shown in Fig. 3.5(a). Similar profile for nitrogen absorption is shown in 
Fig. 3.5(b).Total oxygen and nitrogen pickup, as indicated by the Figs. 3.5(a) and (b), 
is less than 2 ppm and can be ignored, when compared to absorption by other 
mechanisms (as discussed below) 




(b) Magnitude of nitrogen and oxygen absorption through the surface of steel 
bath 

Similar to the procedure followed in the previous section, the increase in oxygen 
content, , during the time step At can be calculated from Eqn (3.11) 

A[6]s=ko^([0]L-[OOAt 

=ko-A([0]L-[Or-)At 

The value of height of the liquid in ladle, Lh, is 600 cm. The values of 
obtained from Eqn (3.1) and (3.7). Hence, 

A[6]s=ko^([0]L-[0]S)At 

AVl 

= 0.0025— (0.0337-0.0715)x0.1 
600 

= 1.51 X 10'*% 

The increase in nitrogen content, A[N]ts , during the time step At can be calculated 
from Eqn (3.12) 

A(N]s =k„.^([N]L-[N]S)At 

A 

For the liquid bath in a ladle, — ^ . Therefore, 

AVh 

A[Nls = k „ ([N];., - [NlS ) At 



1 


0 


(0 f)4ns-n no66)x0.l 


1.4125 X !()■'' % 


(c) Magnitude of nitrogen and oxygen absorption due to dissolution of the 
entrained gas bubbles 

The increase in oxygen content during the small time step At can be calculated from 
Eqn (3.26) 


a[6] 


0.00437tl^_V^ 


bbl 


w. 


At 


And, the increase in nitrogen content during the small time step At can be calculated 
by Eqn (3.25) 

.r.-. 0.0144 t,^^ 

A[N]bb, = ,,, ' - - At 

The different parameters which are needed so as to calculate A[0]i,y and A[N]^y 
be determined from the Eqns (3.13-3.24). Values of parameters used in the calculation 
are given below 



Calculated 


Calculated 

Parameters 

Value 

Parameters 

Value 

At 


m 

1.57 X 10^ 

Uo 

3.87 ■ 

E, 

3.43 X 10*° 

U 

2090.54 

Rc 

74.31 

ao 

14 

Vg 

1.35 X 10^ 

a 

6.59 

X 

55.19 



86.76 


We 

2.93 x 10' 

0^ 

Fr 

10.89 

£ 

Re 

7.58 X 10' 

ho, 

Wl 

4.167x10' 

hN, 


Hence, the increase in oxygen content during the small time step At is 

loV, 


^[0]bbl — 


0.0043711' V, 


At 


W, 

0.00437x1x1.35x10' 


0.1 


4.167x10^ 

1.416 X 10'^% («14.16 ppm) 


The increase in nitrogen content during the small time step At is 


A[N],„ = 


0.0144ti'h Vg 

W, 


At 


0.0144x1x1.35x10' 

4.167x10' 


0.1 


= 4.665 x 10'^ % (w 46.65 ppm) 

A typical profile of the increase in oxygen content caused by dissolution of the air 
bubbles entrained into the bath versus tapping time is shown in Fig. 3.6(a). Similar 
profile for nitrogen absorption is shown in Fig. 3.6(b). [0] and [N] absorption caused 
by dissolution of the air bubbles causes a significant increase in [O] and [N] contents 
in the liquid bath. It can be seen from the Eqn (3.9) that towards the end of tapping the 
total height of the falling stream is lower than that during the preceding time steps due 
to the filling up of the ladle and consequence decrease of Zt. For this reason, the 



entrainment of air bubbles also decreases with time. Figs. 3.6(a) and 3.6(b) show, 
therefore, a continuous decrease in the magnitude of nitrogen and oxygen pickup at 
different stages of tapping. The absorption caused by dissolution of the air bubbles 
entrained into the bath is significant; the total amount of oxygen pickup is 
approximately 13 ppm and that for nitrogen is 43 ppm (assuming 100% absorption 
efiiciency), for 12 minutes of tapping time for a case of 300 ton liquid steel at Bokaro 
steel plant. 

A summary of nitrogen and oxygen pickup by different mechanism is given in Table 
3.3 The total increase in bath oxygen is plotted as a function of time in Fig. 3.7(a) and 
the total increase in bath nitrogen is plotted as a function of time in Fig. 3.7(b). 

It can be seen from Table 3.3 that approximately 90% of oxygen and nitrogen pickup 
is by the dissolution of gases in the entrained bubbles. The effect of tapping rate on 
total absorption is plotted in Fig. 3.8. As shown in the Fig. 3.8 the absorption of 
nitrogen and oxygen increases with increase in tapping rate. Absorption increase may 
be due to the extra turbulence induced in the metal stream and in metal bath owing to 
high tapping rate. There is an optimum tapping rate to balance temperature loss vis-a- 
vis gas pickup. 




distance from nozzle exit, cm distance from nozzle exit, cm 


Fig. 3.5(a): Absorption of [O] Fig. 3.5(b): Absorption of [N] 

through the surface of falling stream through the surface of falling stream 



Fig. 3.6(a): Magnitude of [O] from Fig.3.6(b): Magnitude of [N] from 


entrained air bubbles at different entrained air bubbles at different 


tapping times of liquid steel 


tapping times of liquid steel 



Magnitude of [0] pickup, ppm 



Fig. 3.7(a): Magnitude of [O] pickup 
by all mechanisms at different 
tapping times of liquid steel 



Fig. 3.7(b): Magnitude of [N] pickup 
by all mechanisms at different 
tapping times of liquid steel 



Fig. 3.8: Effect of tapping time on N 2 and O 2 pickup 



Table 3.3: N 2 and O 2 pickup by different mechanisms 



Absorption Mechanisms 

E.xtcnt of 

Absorption 

(ppm) 

Percentage 

of 

Absorption 

O 2 Pickup 

Diffusion through tapping stream 

1.48 

10.17 

Diffusion through the surface of the metal 

pool 

5.0x10"® 

3.45 X 10"^ 

Dissolution of Bubbles 

13.07 

89.83 


N 2 Pickup 

Diffusion through tapping stream 

1.73 

3.86 

Diffusion through the surface of the metal 

pool 

2.0x10"® 

4.46x10"® 

Dissolution of Bubbles ■ 

43.07 

96.14 



3.5.2 Calcination of initial amount of dcoxidant to be added for a given 
composition and temperature of steel 


In this study, our objective is to predict the amount of aluminum to be added to the 
ladle for known initial oxygen dissolved in the steel bath (in ppm). For the case of 300 
ton liquid steel at 1923 K, initially containing 800 ppm dissolved oxygen, Table 3.4 
shows the effect of aluminum addition on the final value of dissolved oxygen ([O]). 
This is for an ideal case when the recovery of aluminum added to the bath is assumed 
to be 100%. 

A user fhendly computer program has been developed to calculate the amount and 
size distribution of alumina inclusions formed by aluminum deoxidation. Flow chart 
of the algorithm has been shown in Fig. 3.9. Various steps of calculation are as 
follows 

Step 1 : Enter relevant input through user fiiendly interface, and save to input file 
“Data.txt” 

Step 2: Set initial time as to zero i.e., t = 0 and the small time step as At = 0.1 sec 
(say). 

Step 3: calculate the equilibrium concentrations of aluminum and oxygen by iteration 
ofEqns (3.34-3.36). 

Step 4: Calculate the super-saturation of the bath by using Eqn (3.46) 

Step 5: Calculate critical radius for nucleation of alumina by using Eqn (3.47) 

Step 6: Calculate bath concentrations of aluminum and oxygen can be calculated by 
Eqn (3.56) and Eqn (3.57) 

Step 7: Calculate interfacial concentration of aluminum, C'^,^, , by using Eqn (3.55) 



Step 8: Caicuialc llie in ladms of uliumna particles by mass transfer of 

aluminum and oxygen to the metal-alumina interface in the small time step At by 
using Eqn (3.59) 

Step 9: Calculate the decrease in bulk concentration of aluminum and oxygen 
(Cai .Co , respectively) due to alumina formation, by using Eqns (3.61) and (3.62) 

Step 10: Calculate the growth and removal of alumina particles by using Eqns (3.44), 
(3.53) and (3.54) 

Step 1 1 : Calculate the individual and total number of particle removed within time 
step At by different collision mechanisms and by assimilation in top slag, by using 
Eqns (3.49), (3.50), (3.53), (3.54) 

Step 12: Store the data for the calculation of time, t, in different output files (see sec 
3.6) 

Step 13: Increment the value of t by At 

Step 14: If C A, > Ca, , then go to ‘Step 4’ else continue 

Step 15: If t < Total time, then go to ‘Step 8’ else continue 

Step 16: Display the results of inclusion size and amount, dissolved oxygen, nitrogen 
and aluminum, and maximum size of inclusion in the output window. 




Fig. 3.9: Flow chart developed for calculation of critical radius of nucleus, growth of 
inclusions, and removal of alumina inclusions, during and after aluminum deoxidation 













Table 3.4: [O],- and final [Al]f vs. Aluminum addition for SOOppm initial [O], 
Recovery 1 00%, temperature 1 873 K 


. No. 

Al added 

final wt 

[0]f 


(kg) 

(ppm) 

(ppm) 

1 

200 

‘1.5 

209.0 

2 

220 

2.5 

150.6 

3 

240 

5 

93.6 

4 

260 

16 

43.9 

5 

280 

55 

19.2 

6 

300 

113 

11.9 

7 

320 

176 

8.8 

8 

340 

241 

7.2 

9 

360 

307 

6.1 

10 

380 

373 

5.4 

11 

400 

439 

4.8 

12 

420 

506 

4.4 

13 

440 

572 

4.0 

14 

460 

639 

3.7 

15 

480 

705 

3.5 

16 

500 

772 

3.3 



'I he final dissolved oxygen in ppm and final mass % aluminum dissolved into the bath 
after deoxidation versus aluminum addition in kg is plotted in Fig. 3.10. As seen from 
the Fig. 3.10, for the case of 800 ppm initial dissolved oxygen, with increasing the 
addition of aluminum the final amount of dissolved oxygen decreases initially at a 
rapid rate, but after 300 kg of aluminum addition (corresponding to equilibrium 
oxygen content of 7 ppm), further decrease in the dissolved oxygen is at a small rate. 
The final oxygen content of the metal bath versus aluminum added into the bath for 
different initial oxygen contents is also plotted in Fig. 3.10 



Fig. 3.10: Plot of A1 used vs. final mass % of dissolved aluminum and oxygen 




3.5.3 Recovery of aluminum during deoxidation 


Percent recovery of aluminum is defined as 


% Recovary = 


Aluminum used for deoxidation + Aluminum dissolved in the metal bath 
Total aluminum charged into the bath 


K 100 


It is clear that aluminum reacting with FeO or O 2 in atmosphere is a loss. 


(a) Calculation of Recovery of aluminum charged during tapping of steel at 
Bokaro steel plant (300 ton steel) 

Heat No: 1442 
Date: 27/2/2004 

Venue: SMS-II, Bokaro Steel Pl ant 


During Tapping: 

Initial Oxygen in the bath [0]i = 506 ppm (measured by Celox) 
C charged = 300 kg 
A1 Charged = 460 kg 


At SRU: 

Final Oxygen in the bath [0]f: < 3 ppm (measured by Celox) 

Final A1 in the bath [Al]f = 0.07 % 

Calculations: 

300 kg C reduces 100 ppm [0] from the bath (approximately 10% reaction assumed). 
Rest 406 ppm [0] was reduced by AI. 



Hence [OJ reduced by A1 


().()4()3x 3x K)^ 

Too 


kg [O] 


=120.9 kg [O] 


Aluminum used = 


^ 2x27 ^ 0,0406x3x10^ 
v3xl6j 100 


= 136.013 kg 


Suppose the initial aluminum content of bath is zero. When the ladle is brought to 
SRU the bath contained 0.07 % dissolved aluminum. 


0 07 x 3 X 1 0^ 

Mass of aluminum dissolved in the metal bath = — ^ — kn fAll 

100 

=210kg[Al] 

Total aluminum used in metal bath = (136.013 + 210) kg 
=346.013 kg 

Since total aluminum added was 460 kg, the percentage recovery of aluminum is 

T> . 346.013 

Percentage recovery = x 100 

460 


= 75.22 % 

Such a high recovery is possible only when slag stopper is used in EOF to prevent 
transfer of slag from EOF to ladle and also coke is used for initial deoxidation 
purposes. 

According to the model, by considering 75% recovery, the equilibrium mass % of 
aluminum is 0.066 % and equilibrium oxygen content in the bath is 3.66 ppm. 
Therefore, it can be concluded that the process has achieved near thermodynamic 


equilibrium. 



(b) Calculation of Recovery of A1 charged at SRU: 

Heat No: 2494 
Date: 10/4/2004 

Venue: SMS-II, Bokaro Steel Plant 


After arrival at SRU: 

Initial oxygen content of bath [O]; = 12.2 ppm (measured by Celox) 
Initial A1 content of bath [Al]i = 0.006 % 

A1 Charged = 180 kg 


After arcing and purging at SRU (Ladle furnace); 

Final Oxygen in the bath [0]f: = 4.3 ppm (measured by Celox) 
Final A1 in the bath [Al]f= 0.016 % 

Calculations; 


[0] reduced by A1 = 


(0.00122 - 0.00043) x 2.893 x IQ- 
100 


kg [0] 


=2.285 kg [O] 


Aluminum used = 


/' 2x27 ^ 
y 3 X 16 j 


X 2.285 kg 


= 2.571 kg 

Again, in the bath %A1 increases to 0.016% from 0.006% 


„ (0.016 -0.006)x3x 10^ , 

Hence A1 dissolved into the bath kg [Al] 


= 30 kg [Al] 


Total Al used in metal bath = (2.285+30) kg - 32.285 kg 



n. 32.285 

Hence, % recovery = xlOO = 17.94 % 

1 80 

The recovery of aluminum at ladle furnace is small because a large proportion of 
aluminum, added in the form of aluminum shots, melts at surface and reacts with 
ambient atmosphere. Argon stirring also helps in oxygen pickup from atmosphere. 
The aluminum wire injector is not used at Bokaro steel plant. If aluminum wire is 
practiced then recovery will improve. 

According to the model, by considering 18% recovery, the equilibrium mass % of 
aluminum is 0.0165 % and equilibrium oxygen content in the bath is 9.26 ppm. 


3.5.4 Calculation of initial size of nucleus and its growth by mass transfer 

The initial conditions assumed for calculating the size of the nucleus (from Eqn 3.55) 
and its subsequent growth (from Eqns 3.61-3.65) are summarized in Table. 3.5. 

(a) Calculation of super-saturation value and critical radius 


The super satxoration, 11, can be calculated from Eqn (3.45) 



The values of Ce and C, can be obtained by calculating equilibrium constant for the 
deoxidation reaction from Eqn (3.33) 


log(Ke) = 


58473 

TK 


-17.74 


17.74 

1923 


= 12.66718 


Table 3.5; Data used for the calculation 


Symbol 

Meaning 

Unit 

Numerical 

Value used 

Ali 

Initial A1 addition into the bath 

i 

kg 

400 

dpe 

Average density of iron ~ 

— — — 

kg/m"* 

7000 

din 

Average density of inclusions 

kg/m'^ 

3000 

H 

Height of liquid bath of the ladle 

m 

3.5 

Kai 

Average difflisivity of [Al] in molten steel 

m'^/s 

2.5 X 10’^ 

Ko 

Average diffusivity of [0] in molten steel 

m'^/s 

2.5 X 10'^ 

No 

Initial no of particles 


1.0 X 10'‘^ 

[0]i 

Initial Oxygen level 

ppm 

800 

R 

Universal gas constant 

J/K 

8.314 

tn 

Total number of moles present in the bath 

mole 

5357.14 

T 

Temperature of steel in Kelvin 

K 

1923 

Vn, 

Volume of metal 



42.86 

Wn, 

Weight of metal 

kg 

3 X 10" 

TlAl 

Aluminum recovery 

% 

100 

|iFe 

Average viscosity of iron 

kg/m.s 

0.005 

CJs-i 

Surface tension at steel and inclusion interface 

■■ 

2.3 






Or, PCe = = 4.647U4 x 10‘‘ 

By using Eqn (3.34 - 3.36), the values of [Al]e and [0]e are 0.04395 % and 0.0004S 
%, respectively. Hence 

Ce=[Al]e^[0]e^ = 2.15x 10-'^ 

The values of [Al]i and [0]i are 0.13333 % and 0.07 %, respectively, for the first time 
step. Thus, 

Ct = [Al]i^ [0]i^ = 9.102223 x 10'^ 

The super-saturation of the bath is, 
log (Ct/Ce) = 7.62639 

The critical radius of nucleus can be calculated from Eqn (3.47) 

_ 2CTVm 

X- 

RTlog ^ 

V^e 7 

^ 2x2.3x3.4x10-^ 

8.314x1923x7.626391 

= 0.1283 X 10'^ m 

0 

Critical radius is approximately 0. 1 3 A 

The calculated value of the critical radius calculated from classical homogeneous 
nucleation theory is less than the size of atom. This is primaiily due to large super 
saturation assumed. For further calculation the value of critical radius is assumed to 

o 

be 8 A , which would be approximately the value when Ct is only slightly greater than 


Ce (for example log 



0 . 1 ) 


(b) Growth of nucleus by mass transfer 


When the deoxidation achieves equilibrium then the for a small time step, At the 
increase in radius, Ar, can be calculated from Eqn (3.59) 

Ar= 0.01275 Rai (C^, -Ci,,) At 

Here, the values of C^j and C^| can be calculated from Eqns (3.55) and (3.56). 

From Eqn (3.55), 


=[mass%Al’’] — ^ — 

100 xM^, 


Where, Mai is the molecular weight of aluminum, and is density of iron. 

7000 

^100x27 

= 0. 3457 kmole/m^ 


Hence, C*, = 


400 

300000 


-xlOO 


The bulk concentration of oxygen can be calculated from Eqn (3.57) 


C' =[mass%0^] — — 

° lOOxMo 

Where, Mo is the molecular weight of oxygen and is equal to 16. 

Hence, Ct = 0.07 x = 0.30626 kmoleW 

° 100x16 

The value of C‘^, has been calculated by iteration method from Eqn (3.55), the value 
of C’,^, for the first time step is 0.1429 Ehnoles/m^. 

Now, the value of interfacial concentration of oxygen can be calculated from Eqn 
(3.50) 


1 


C- 


f 


"Fe 




lUU.-.M,, ) 


K. 


C'a, X 


Pai 

VP.X 


^ 7000 ^ 

^100x16; 


1 

i 4.64704x10'^ X 

'0.142942200^1 

1 

L V7000 JJ 


2 


= 0.00181 kmole/m^ 


Thus, increase in the radius during a small time step At is 
Ar = 0.01275 kA!(C^, - Q,)At 

= 0.01275 X 2.5 X 10'^ (0. 3457 x 10'^-0.1429 x lO'^) x 0.1 
= 6.464 xlO'^^ m (or, 6.464 xlO'^ pm) 

Corresponding increase in mass of alumina in the bath after deoxidation can be 
calculated from Eqn (3.60) 


d[AL03] = No— Tc{[r + dr]^ r }Pai203 

Where, No is the initial number of inclusions formed in the bath and assumed to be 
equal to lO'^. Now by putting in the values of No, r ,dr and Pai203 > 


d[Al2O3] = 10’^ x^7r{[8xl0'* + 6.464 xl0'’^f -[8x10 3000 

= 1.56x10'® kg 

Change in mass percent of aluminum in bulk can be calculated from the Eqn (3.61) 


„ 2xd[Al203] 

d[mass% Al] = ,, ' — x 1 00 

W, 

- 2x1-56x10 .^|nno/. = 1.04 x 10''^% 

. 300000 

Change in mass percent of oxygen in bulk can be calculated from the Eqn (3.62) 




3x1.56x10’ 

= -xlOO 

300000 

=1.56 X 10'’^% 

The radius of inclusion is plotted as a function of time in Fig. 3.11(a) and the 
concentration of dissolved oxygen and aluminum are plotted as a function of time in 
Fig. 3.11(b). It is clear from Fig. 3.11(a) and 3.11(b) that the system attains chemical 
equilibrium after approximately 5 minutes of aluminum addition and the minimum radius of 
inclusions present at equilibrium is approximately 2.4 microns. The final size would increase 
if number of nuclei is assumed to be less. 



time, min time, min 

Fig 3. 1 1 (a): Decrease of mass % [Al] Fig 3.1 1 (b): Increase of inclusion radius 

and [O] before equilibrium before the system achieves equilibrium 





3.5.5 Calculation of change in the size range and number of inclusions by 
collision and assimilation by top slag 


The equations to be used for calculating the change in size and number of inclusions 
due to collision and assimilation by top slag have been described in section 3.4. The 
simplifying assumptions made in the calculation, procedure adopted in this work, are 

- Deoxidizers dissolved uniformly in the melt soon after addition. 

- The inclusions formed are also homogeneously distributed through out the steel 
melt. 

- A certain initial number of nuclei are assumed to be formed in the bath 

- After attainment of thermodynamic equilibrium, the inclusions grow only by 
collision mechanisms and not by mass transfer because system has already achieved 
equilibrium 

Table 3.5 summarizes the data used in the kinetic model for the prediction. The 
sample calculations are given below. 



Table 3.5: Data used in the kinetic model 


S>Tnbol 

Meaning n 

Unit 

Numerical 




Value used 

Ali 

Initial A1 addition into the bath 

kg 

400 

dpe 

Average density of iron 

kg/m"* 

7000 

din 

Average density of inclusions 

kg/m"^ 

3000 

H 

Bath Height ~ ^ 

m 

3.5 

K 

Boltzmann constant 

J/K 

1.38 10'^^ 

Kai 

Average diffusivity of [Al] in molten steel 

m"/s 

2.5 X 10'^ 

Ko 

Average diffusivity of [0] in molten steel 

m^/s 


No 

Initial no of particle 


1.0 X lO'"^ 

[0]i 

Initial Oxygen level 

ppm 

800 

R 

Universal gas constant 

J/K 

8.314 

T 

Temperature of steel in centigrade 

°C 

1650 

TK 

Temperature of steel in Kelvin 

K 

1923 

Vm 

Volume of metal 

m'^/raol 

3.4 X 10'^ 

Wm 

Weight of metal 

kg 

3x 10^ 

a 

Collision efficiency 


0.5 

riAi 

Aluminum recovery 

% 

100 

M’Fe 

Average viscosity of iron 

kg/m.s 

0.005 

1 

1 

Surface tension at steel and inclusion interface 

N/m- 

2.3 


The calculation is done over a small time step. At = 0.1 sec. The collision volumes of 
different mechanisms can be calculated from the Eqns (3.63-3.65) A sample 
calculation for the first time step is given below. 

From the calculation of equilibrium radius of inclusions, the minimum size of 
inclusion present in the bath after attainment of equilibrium, reqim, is approximately 
equals to 2.4 pm. The time required for attainment of equilibrium, is 

approximately 4 minutes and the total number of inclusions present in the bath, No ,is 
assumed to lO’^. 


at t = 0, 

Ni = No=10'^ifi=l 
Ni = 0ifi = 2, 3, 4 


Now, at t = 0 + At 

The collision volume due to Brownian motion can be calculated from Eqn (3.63) 


3p 


Fe 






2x1.38x10“"^ X1923N 1 


(-+-)(^o+^o) 


^0 h 


3x0.005 

= 21-3810-^^x1923 
3x0.005 

= 3.53832x 10-’* 

The collision volume due to Stokes collision can be calculated from Eqn (3.64) 




^^{pFe -Pi) 


V 






Fe 


As here, i = j =1 [Only one size range particle exists at t - 0] 
Hence, r,- = rj = ri = ro 



w%, = 0 

So, no Stokes collision takes place during the first time step. 

Calculation of decrease in the number of inclusions by collision mechanisms, as 
described by Eqn (3.67), is 

dNi (i == 1) = - (Wii^ + Wii®) X Ni xNi X dt 
= - [0 + 3.53832 X lO’’* ] x 10^^ x lO’^ x 1 
= - 3.53832 X 10'^* X 10'^ X 10'^ 


= -3538320 

dNj C = 1) = - (Wii^ + Wii®) X Ni X Ni X dt 
= -3538320 

dN[i+i] = Na = + (Wi + Wi 1®) X N, X N, X dt 
= 3538320 

Therefore, Ni*= Ni + 2 x dNi = 10*^- 2 x 3538320 = 99999292340 
N2*= Na + 2 x dNa = 0 + 2 X 3538320 = 7076640 
Hence for the first step the decrease in number of inclusions due to collision between 
particles is = dNi +dNi -dNa 


= 2x3538320-3538320 
= 3538320 


Now, in the same time step some of the ra particles will collide with both ri and ra 
sized particles and they will form inclusions of size ra and r 4 , respectively, and so on. 
The decrease in number of inclusions removed by assimilation by top slag can be 
calculated from Eqn (3.68) 

dt 9//^^^^ H 


For i =1, 



^ 2g 

di 9/; 



// 


Or, 

9/J. H 


2x9.81 

9x0.005 


(7000 - 3000)(2.37772088 x 1 0'^)^ 


999962509312 , 

xl 

3.5 


= 281698 


Similarly dN 2 = 8 , dNs = 0, dN 4 = 0 and so on, for the first time step. 

Hence the total removal to the slag = dNi + dNa + dNs + dN 4 

= 281698 + 8 + 0 + 0.... 

= 281706 

Result obtained after a time step, At, is 
Ni = 999962509312 - 281698 = 999962227614 
N 2 = 1 8736086 - 8 = 1 8736078 
N 3 = 0 , N 4 = 0 and so on. 

These values are to be used in the calculation during the next time step. At. 

The logarithmic distribution of inclusions of different sizes, which are present in the 
liquid bath after different ladle holding times, is plotted in Fig. 3.12. As seen firom the 
Fig. 3.12, the range of inclusion sizes present in the bath widens with increase in the 
holding time. The concentration of small sized inclusions in the bath decreases with 
holding time, and the concentration of bigger size inclusions increases at the expense 
of the smaller inclusions. The total number of inclusions continuously decreases 


(Table 3.6 and Fig. 3.13). 
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Fig. 3.12: Plot of logarithmic inclusion distribution for different holding times 
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Fig. 3.13: Predicted number of inclusions and mass of inclusions present after 
different holding times 






l ablc 3.6; Inclusion size distribution in the bath after different holding times 


Range(fim) 

5 min 

10 min 

20 min 

30 min 

45 min 

60 min 

2 -3 

9.00 X 10” 

9.92 X 10” 

9.85 X 10" 

9.72 X 10” 

9.514 X 10” 

9.087 X 10” 

3 -4 

2.33 X 10® 

6.56 X 10® 

1.61 X 10® 

3.04 X 10® 

4.5 X 10® 

5508048384 

4 -5 

4788148 

39889460 

2.5x10® 

9.52 X 10® 

2.167 X 10® 

3597425664 

5 -6 

981 

111649 

5709985 

86377520 

436471136 

1206495616 

6 -7 


87 

76555 

6673791 

89937752 

444622976 

7 -8 



630 

452026 

19684498 

190102880 

8 -9 




16176 

3153790 

68448688 

9 -10 




472 

497124 

26183152 

10 -11 




10 

70916 

10014960 

11 -12 





8706 

3607955 

12 -13 





1083 

1385015 

13 -14 





121 

515376 

14 -15 






191154 

15 -16 






71870 

16 -17 






28210 

17 -18 






10941 

18 -19 






4434 

19 -20 






1800 

20 -21 






755 

21 -22 






335 

22 -23 






149 

23 -24 






70 

24 -25 






34 

25 -26 






19 

Total 

9.96 X 10” 

9.93 X 10” 

9.87x10” 

9.76 X 10” 

9.59x10” 

9.20x10” 



The total number of inclusions and the mass of inclusions present in the bath is shown 
for different holding times in Fig. 3.13. It is assumed here that turbulent stirring is 
absent during the time of transportation of the ladle to SRU. According to Fig. 3.13 in 
absence of stirring only 2 % of the inclusion formed by the deoxidation reaction are 
likely to be removed within first 20 minutes of ladle holding time, and 7% in 45 
minutes of ladle holding time. For higher ladle holding times the chances of formation 
of bigger particles increase (Fig. 3.14) due to the contribution of stokes collision. 
Also, the chances of removal of inclusions by flotation increase as the size increases. 
The largest size inclusion formed in the bath is plotted as a function of time in Fig. 
3.14; after approximately 5 minutes of holding time, there is a change in slope of the 
curve. This is due to the assumption that the contribution of Stokes’ collision begins 
only after the inclusions attains a size of 5pm. It has been reported that for the case of 
Stokes collision the collision volume decreases rapidly in the size range of 2-4 
microns Hence, during the first 5-6 minutes Brownian motion is, perhaps, the 
dominant growth mechanism. After 5-6 minutes, there is a sufficient number of 
inclusions in size range of 4-6 microns, and therefore the Stokes collision mechanism 
becomes the dominant growth mechanism. This is also evident from Fig. 3.15. 
According to the model, it is possible to form inclusions of size of 15 pm in the metal 
bath by holding but without purging. It can be seen fi'om Table 3.6, the number of 
bigger inclusions is less, but the number of inclusions of size 5-10 pm is high. A 
sample of Bokaro steel plant (BSP), collected before purging (approximately 45 
minute holding time) may confirm this. A SEM study of sample collected at SRU just 
before purging is shown in Fig. 3.16. The diameter of a typical inclusion, in Fig. 3.16, 
is approximately 15 pm. 




Fig. 3.14: Prediction of the largest sized inclusion in the bath for different holding 
times 



Fig. 3.15: Contribution of different mechanisms to the total removal of inclusion 






Fig. 3.16: SEM micrograph of a sample (before pmging) collected from BSP showing 


a typical alumina inclusion of approximately 15 pm size. 



time, min 

Fig. 3.17: Rate of change of different inclusion removal mechanisms with time 



Fig. 3.15 compares the contribution of different mechanisms to inclusion removal. It 
can be seen that in the beginning, the contribution of Brownian motion dominates, but 
later on Stokes collision becomes a more dominant mechanism of inclusion removal. 
The rate of inclusion removal by the collision caused by Brownian motion is nearly 
constant (infact, it decreases at a slow rate within the 45 minutes of holding time), as 
shown in Fig. 3.17. It may be due to the fact that the decrease of smallest sized 
inclusion is only 8.3% of the initial concentration within a 45 minute holding time. 
The rate of assimilation of inclusion by the top slag slowly increases with time. After 
5 minutes of holding time the rate of inclusion removal rate by Stokes’ collision 
increases rapidly fi-om a very low value. 

3.5.6 Effect of gas purging on the total removal of inclusions 

A study of the relative contribution of different collision mechanisms on total removal 
of alumina inclusions in the presence of gas purging is presented in this section. The 
comparison of contribution of mass transfer and other mechanisms to inclusion 
growth, before the bath attains equilibrium, is shown in Fig. 3.18. It can be seen fi'om 
Fig. 3.18 that the largest size inclusion formed in the bath by mass transfer 
mechanism gradually increases with time, whereas, by introducing other collision 
mechanisms Rmex increases rapidly (with sudden jumps). The reason for sodden 
jumps is due to the formation of bigger particles as a result of collision at discrete 
time steps, as assumed in the calculation procedure in this work. 




Fig. 3.18: Comparison of increase in Rwax when mass transfer mechanism and when 
all mechanisms are acting in the first 4 minutes of deoxidation 



Fig. 3.19: Comparison of increase in Rmx due to only mass transfer, due to mass 
transfer and collisions (without purging), and due to all mechanisms act before the 
bath attains equilibrium 




It can be seen from Fig. 3.19 that the effect of purging (of 2 lit/ton nun gas) on Rmox is 
only limited. But the presence of collision mechanisms before equilibrium affects the 
value of RMax significantly. This may be due to the contribution of turbulence in the 
bath by purging of gas. Stirring enhances the chance of collision and hence the chance 
of formation of bigger particle. Therefore, the contribution of other mechanisms also 
increases (Figs. 3.20(a) and (b)). 

The effect of gas purging on the rate of removal by stokes and Brownian motion is 
shown in Fig. 3.20(a) and (b). It can be seen from Fig. 3.20(a), that due to 
introduction of turbulence into the bath (2 lit/min/ton of metal), the rate of inclusion 
removal by stokes collision increases by a factor of 10 compared to the case when no 
purging is used. It can also be seen from Fig. 3.20(b), that the effect of purging on the 
removal rate of inclusions by Brownian motion in not significant, which is expected. 
The comparison of total inclusion removal from the bath, with or without purging is 
plotted in Fig. 3.21. Purging with a gas flow rate, Q, of 2 lit/min/ton of metal is 
capable of reducing inclusions by approximately 50% (of the initial number of 
inclusions) in less than an hour of holding time. 



% of inclusion removed dNs/rnin 



Fig. 3.20(a): Effect of purging 
on Stokes collision 



time, min 

Fig. 3.20(b): Effect of purging 
on Brownian motion 



Fig. 3.21 : The effect of turbulence collision on total removal of inclusions 




3.6 User friendly model for prediction of size distribution and 
amount of alumina inclusion for aluminum killed steels 


A prediction model, incorporating the kinetics of aluminum deoxidation, various 
collision mechanisms and assimilation by the top slag, has been developed. The 
sequence of steps executed in the model is shown in Fig. 3.22 




Saved input 
Data 


Control 

Model 

User interface 







Activated 





Display 



Save 

Results 



Results 


Fig. 3.22: The sequence of steps executed in the user friendly model 


It may be noted that in the initial stages of deoxidation (say during the first four 
minutes) the inclusion size range is small, owing to assumption made that all nuclei 
are of the same size to start with. Due to this reason the time needed for computation 
of inclusion growth and coalescence by different mechanisms is also small. The 
computation time, however, increases almost exponentially for large holding times. A 
special optimization procedure based on genetic algorithm (described in Appendix 
3. A) has been incorporated in the present work, for the first time, to reduce the 










computation time, for example by nearly 70%, for holding time of 60 minutes or 
more. 

The user friendly interface of this prediction model is shown in Fig.3.24. The features 
of the user friendly interface are as follows: 

• Input window: This facility helps the user to key-in the initial conditions, 
including the dissolved oxygen measured in the ladle at the time of sampling, 
estimated nitrogen content, amount of aluminum added for deoxidation, 
tapping time, holding time, temperature, etc. 

• Output window: Displays the results of inclusion size and amount, dissolved 
oxygen, nitrogen and aluminum, and maximum size of inclusion. 

• Storage facility: Stores data in input and output folders which can later be 
retrieved for further analysis and verification. 

• Graph: Model provides online graph facility. The plot of inclusion 
distribution immediately after the calculation can be obtained. 

3.6.1 Steps to be followed for Using the Model: 

Stepl: Double click on the desktop icon ‘ Collision, exe’. 

Step 2: The “input window” will open. The inputs can be changed by typing-in a new 
value in the corresponding text box. Click ‘ENTER’ and the out put window 
(Fig.3.25) will appear. 



Step 3: For viewing the plot of inclusion size distribution click on the “PLOT” 
button. It will generate the plot of size range vs. logarithm of number of inclusions 
(see Fig. 3.26). 

Step 4: For details of results and further access see following output files created by 
the program (inside the installed folder). The output of different output files are as 
follows 

• Absorption.txt: Gives the predicted amount of oxygen and nitrogen absorbed 
by the liquid steel during tapping 

• B4EqIm.txt: Gives the increase in inclusion radius, decrease in dissolved 
aluminiun and oxygen content before attainment of equilibrium in the bath 

• SizeCaL.txt: Gives the predicted amount of inclusion of different sizes after 
total holding time 

• MaxR.txt: Gives predicted maximum radius for different times 

• Removed.txt: Gives the predicted amoxmt of inclusion removal by different 
mechanisms 

• Outputtxt: Gives the predicted amount of inclusion of different size ranges 
after total holding time 
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Fig. 3.24: User friendly input window of the model 
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Fig. 3.23: Desktop icon of the user jBnendly model 
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Fig. 3.25; Output window of the model 
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Fig. 3.26; The plot of amount of inclusion vs. range of inclusions in the output 


window of the model 





4. Conclusions 


4.1 Conclusions for chapter 2 

A good strategy, after having obtained the operational data, is to do statistical analysis 
of the data in terms of normal distribution, clustering, covariance, etc. and then do 
multiple linear regression analysis to develop regression models. Simultaneously, 
kinetic model based on first fundamentals should be developed so as to understand the 
process well. This will eventually help in deciding the direction of any technical 
improvement/modification. The next step should be to optimize the kinetic model by 
GA, develop ANN model as a back-up support, and finally, if needed, develop 
FREES model for predictions to be made only in very specific process operation 
ranges. The results of different models, when analyzed together, help in cross 
checking the performance of models in different regions of operation. In a plant 
environment, one should not base all decisions on the basis of one model only but 

rather narrow down to the possible decision range and then take an optimal decision 
[ 2 . 22 ] 

4.2 Conclusions for chapter 3 

For 300 ton liquid steel, having initial dissolved oxygen of 800 ppm and nitrogen of 
50 ppm, and deoxidized with 400 kg aluminum 

1. Absorption of [O] and [N] through the surface of the falling stream may 
contribute approximately 7-8 % of the total pickup of [O] and [N]. The entrained 
gas bubbles contribute approximately to 90 % of the total pickup of [O] and {N]. 



2. In order to obtain a final oxygen content of bath below 5 ppm the optimum 
aluminum to be added lies between 350 to 500 kg, depending on percentage 
recovery of aluminum addition. 

3. In the case of non-stirred bath, the concentration of inclusion decreases with time 
at a low rate. A 45 minute holding time is capable of removing only 7% of the 
initial inclusions present. Approximately 10% of the initial concentration of the 
smallest inclusion, Rmin, can be removed within 45 minute holding time. 

4. The contribution of Stokes’ collision rapidly increases with time as the 
availability of the big inclusions also increases with time. Similarly, the 
contribution of assimilation of inclusion by top slag also increases with time 
because the big inclusions have high upward velocity. 

5. During holding of steel in a non stirred-bath, the collision mechanism plays the 
most significant role for removal of inclusions. For a 45 minute holding period 
the inclusion removal by collision mechanism is more than 15 times than that by 
assimilation of inclusions by slag. 

6. In case the bath is stirred by gas or induction current turbulent, collision plays a 
significant role on the inclusion removal process. For example, if the bath 
isstirred with a gas flow rate of 2 lit/ton/min, approximately 50% of the initial 
inclusions may be removed in 45 minutes. 

7. Inclusion removal process is a dynamic process hence the time step for 
calculation should be small (0.15). Results obtained by using high At value may 


lead to higher error level. 
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Appendix-A 


Optimization of computation time by genetic algorithm 

The usability a model in an industry is depends on the computation time required for 
execution of the model. It has been foimd, in this study, that at the beginning (say for 
time less than 20 minutes) the total number size of inclusions are less, and for higher 
holding time (say more than 45 minutes) that is large. The size range and the 
corresponding number of inclusions present in the bath for different holding times, as 
predicted by the model, are presented in Table 3.7. 

Table 3.7: Number of size range present at different time 
No. of sizes 

Time, min (n) 


2 

8 

5 

12 

10 

19 

20 

37 

30 

70 

45 

200 

60 

751 


If a large size range is assumed from the beginning of the calculation, the computation 
is extremely slow, even for 20 minutes holding time. To reduce the computation time 
linear regression and genetic algorithm has been applied separately. A linear 
regression model has very high value (R^= 0.99). It can be seen from Fig. 3.27, 
linear regression gives a closer result to the data shown in Table 3.7. The linear 
regression model, for some time step, predicting lesser size range than that of the 
model (for example 2"'^ minute, 30'*’ minute) a.s shown in the Fig. 3.2 ■ 1 herefore. it 



introduces compulation^! errors a! different time steps as the predictions may lesser 
than that of the model. Optimization by genetic algorithm (G.A.) gives better result. 
The optimization technique with respect to the data, given at the Table 3.7, is 
described below 

It is assumed that the size range is follows a 3'^'* degree polynomial function 
F(x) = A t^ + B t^ + C t +D. 

Where, A, B, C, D are the constants and needed to be find out by G.A, t is time, and n 
is the number of sizes at different time steps as calculated by the model. Predicted 
F(x) should be maintain the tolerance with the data given by the model, n, 
i.e.; F(x) > n 

For providing sufficient room over n, the constraint has been redefined as 
F(x)>n+ Vst 

Therefore, the single objective optimization problem is as follows 
Minimize F(x) = A t^ + B t^ + C t +D 

Subjected to: F(x) > n + Vs t 

Optimization of the function F(x) by genetic algorithm has been done by using a GA 
code available at the website of Kanpur Genetic Algorithm Laboratory (KAJSIGAL). 
The different GA parameters are given below 

Number of generations is 5000, population size is 300, cross over probability is 0.8, 
mutation probability is 0.2, distribution index for SBX is 2, and distribution index for 
mutation is 100 

The values of the constants as given by GA are A = 0.01, B = - 0.57, C = 12.25, D = 
5.3 1 7. The prediction of number of size ranges of inclusions is shown in Fig. 3.27. 




Fig 3.27; Comparison of predictions by both linear regression and genetic algorithm 
for the prediction of number of size range of inclusions 


The comparison of predictions by both linear regression and genetic algorithm is 
shown in Fig. 3.27. From Fig. 3.27, it can be seen that by using genetic algorithm the 
number of size range is always above the required number of size range, predicted by 
the model, at different time steps. For the first 20 minutes of holding time, the number 
of size range is small, so if GA predicts 50 more size ranges then it will not affect the 
, computation time much. Linear regression, sometimes, predicts lesser number of size 
ranges, and as discussed earlier that introduces calculation error. Thus, it can be 
concluded that by introducing some tolerance, the chances of error is reduced. 




